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PREFACE 


This  report  examines  the  feasibility  of  accounting  for  the  processes  of 
erosion  and  overwash  in  the  determination,  by  the  Federal  Emergency  Manage¬ 
ment  Agency  (FEMA),  of  coastal  high-hazard  flood  zones. 
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APPENDIX  A:  DISCUSSION  AND  ANNOTATED  BIBLIOGRAPHY  OF  OVERWASH 
LITERATURE  . 


Introduction  . 

Selected  Bibliography  on  Overwash 
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PART  I:  INTRODUCTION 


Background 


1.  In  1985,  an  unusual  sequence  of  three  hurricanes  entered  the  Gulf 
of  Mexico  and  moved  toward  Grand  Isle,  Louisiana.  Although  the  newly 
constructed  protective  beach  and  dune  survived  the  passage  of  Hurricanes 
Danny  and  Elena,  Hurricane  Juan,  which  followed,  caused  major  damage  to 
property  and  leveled  large  sections  of  the  remaining  dune.  In  the  same  year 
Hurricane  Gloria  posed  a  major  threat  to  the  east  coast  as  the  first  major 
hurricane  to  pass  the  heavily  developed  mid-Atlantic  shore  since  Hurricane 
Hazel  in  1954.  In  addition  to  hurricanes,  more  frequently  occurring  winter 
storms  regularly  cause  erosion,  flooding,  and  property  damage  along  the 
Atlantic  and  Pacific  Ocean  coasts  (Figure  1). 


Figure  1.  Storm  damage  caused  by  erosion  of  the  dune 


2.  Because  of  the  severity  of  coastal  storms  and  their  potential  for 


causing  loss  of  life  and  property  damage,  understanding  and  prediction  of 
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all  aspects  of  storms  and  their  impact  has  been  a  fundamental  objective  of 

/ 

coastal  research.  Unfortunately,  the  understanding  of  sediment  transport  and 
beach  erosion  lags  behind  the  capability  to  predict  the  waves  and  flooding 
which  may  occur  during  a  storm  event. 

3.  The  National  Flood  Insurance  Program,  which  insures  property 
located  within  coastal  and  inland  flood-prone  areas,  is  based  on  the  expected 
flood  levels  associated  with  a  storm  with  a  1  percent  probability  of 
occurring  per  year  (a  100-year  storm).  Wave  heights  are  considered,  but 
there  is  no  official  procedure  to  treat  flooding  resulting  from  either 
long-term  erosion  or  from  the  erosion  associated  with  the  100-year  storm. 
Although  this  may  be  a  reasonable  assumption  for  riverine  flooding,  it  is  not 
appropriate  along  the  coast  where  the  shoreline  is  constantly  evolving. 

Purpose 

4.  In  December  1985,  the  Federal  Emergency  Management  Agency  (FEMA) , 
which  manages  the  National  Flood  Insurance  Program,  requested  that  the  US 
Army  Engineer  Waterways  Experiment  Station's  (WES's)  Coastal  Engineering 
Research  Center  (CERC)  conduct  a  6-month  study  to  examine  the  feasibility  of 
accounting  for  the  processes  of  erosion  and  overwash  in  determining  coastal 
high-hazard  flood  zones  along  the  Atlantic  and  Gulf  coasts.  Specifically, 
this  report  addresses  methods  for  determining  variations  in  the  width  of  the 
velocity  (V-)  zone,  as  defined  by  FEMA,  which  would  occur  if  these  processes 
could  be  quantitatively  estimated. 

5.  At  present,  an  area  is  designated  a  V-zone  only  if  the  maximum 
water  level  associated  with  the  100-year  storm  would  support  a  3  ft-high 
(0.9  m)  wave.  In  areas  where  dunes  exist,  if  the  dune  height  exceeds  the 
water  level,  the  dune  and  inland  zones  are  assumed  to  be  "protected,"  and  the 
V-zone  is  topographically  mapped  only  to  the  intersection  of  the  water  level 
with  the  dune.  That  the  dune  may  erode  or  that  an  area  is  historically  prone 
to  flooding  during  storms  is  considered  only  by  subjective  means. 

6.  The  fundamental  question  to  be  addressed  is,  "Will  a  particular 
dune  survive  the  100-year  event?"  Specific  details  as  to  how  the  dune 
survives  (or  erodes)  are  not  required.  At  one  extreme  (the  present  FEMA 
procedure),  the  dunes  can  be  treated  as  fixed  barriers,  resulting  in  narrow 
V-zones;  at  the  other  extreme,  the  existence  of  dunes  can  be  ignored, 


resulting  in  wide  V-zones.  The  latter  approach  would  be  easy  to  apply  and  in 
some  areas  may  be  the  only  rational  approach.  In  actuality,  dunes  can 
provide  an  important  level  of  protection  if  they  are  of  sufficient  height  and 
width.  In  those  areas  with  effective  dunes,  a  method  for  estimating 
potential  storm  erosion  is  desirable. 

7.  To  be  of  use  to  FEMA,  a  dune  erosion  prediction  method  must  fulfill 
a  number  of  criteria.  It  must 

a.  Be  compatible  with  existing  FEMA  methodologies  for  computing 
the  nearshore  hydrodynamics. 

b.  Incorporate  FEMA  predictions  of  water  level  and  wave  height. 

c_.  Be  easy  to  use  and  allow  consistent  application. 

d.  Be  relatively  insensitive  to  uncertainty  in  input  variables 
such  as  the  prestorm  beach  configuration,  the  duration  of  the 
storm,  and  the  sediment  size. 

e.  Be  deterministic  and  incorporate  relationships  which  link  the 
erosion  process  to  the  physical  characteristics  of  the  storm 
event . 

f .  Have  undergone  calibration  and  verification  with  field  data. 

Considerable  research  is  currently  under  way  which  will  improve  both  our 
understanding  of  erosion/overwash  processes  and  our  ability  to  predict  them. 
This  study  is  based  on  the  capabilities  of  existing  methodologies.  The  three 
candidate  methodologies  finally  selected  for  detailed  consideration  have  only 
recently  become  available,  dating  from  1982  or  later,  and  work  on  them 
continues.  New  reports  on  each  were  published  subsequent  to  the  completion 
of  the  first  draft  of  this  report  and  are  briefly  noted  here. 


8.  The  report  herein  is  organized  into  seven  parts  and  one  Appendix. 
Part  II  introduces  basic  nomenclature,  discusses  the  process  of  beach 
erosion,  and  introduces  the  concept  of  an  equilibrium  profile.  Part  III 
examines  in  detail  the  requirements  and  suitability  of  theoretical  sediment 
transport  models.  Part  IV  discusses  three  models  which  are  currently  used 
for  estimating  dune  erosion.  Two  of  the  methods  are  compared  in  Part  V. 


Model  selection  and  use  in  V-zone  determinations  are  discussed  in  Part  VI. 
Part  VII  summarizes  the  report  and  makes  recommendations  for  future  research. 
A  review  of  overwash  literature  is  included  in  Appendix  A. 


FART  II:  PHYSICAL  ASPECTS  OF  BEACH  PROFILE  CHANGE 


9.  This  section  discusses  the  general  physical  aspects  of  beach 
profile  change  and  dune  erosion.  The  nomenclature  associated  with  beaches 
will  be  introduced  together  with  the  major  variables  believed  to  control 
profile  change  and  dune  erosion.  The  natural  variability  found  in  measured 
profile  changes  is  quantified,  and  the  equilibrium  beach  profile,  a  funda¬ 
mental  concept  in  present  dune  erosion  models,  is  introduced.  Finally, 
overwash,  a  major  process  associated  with  storms,  is  discussed.  Supple¬ 
mentary  discussions  on  the  topics  of  beach  profile  change  and  littoral 
processes  can  be  found  in  Volume  I  of  the  Shore  Protection  Manual  (SPM) 
(1984) . 


Nomenclature 


10.  Standard  nomenclature  for  describing  beach  profiles  is  shown  in 
Figure  2.  The  locations  and  boundaries  of  various  regions  are  somewhat 
ambiguous  because  of  the  dynamic  nature  of  the  beach.  Definitions  and  names 
of  certain  features  and  regions  may  also  vary  from  author  to  author. 


Figure  2.  Profile  features  and  terminology  (SPM  1984) 


11.  The  profile  is  usually  measured  at  intervals  from  some  point 
landward  of  the  dunes,  across  the  profile,  until  a  depth  is  reached  where 
little  net  sand  movement  is  expected.  This  seaward  limit,  or  "closure 
depth,"  is  on  the  order  of  6  to  8  m  (Hallermeier  1979,  Birkemeier  1985b). 

For  convenience,  the  beach  profile  is  conceptually  divided  into  a  subaqueous 


zone  and  a  subaerial  zone.  The  two  zones  shift  with  changes  in  water  level. 

12.  The  subaqueous  portion  of  the  beach  profile  is  in  continuous 
motion  under  hydraulic  and  mechanical  forces  produced  by  waves  and  currents. 
The  regions  of  most  intense  sediment  movement  are  located  in  the  surf  zone, 
from  the  location  of  wave  breaking  to  the  point  of  maximum  wave  uprush. 
On-offshore  sediment  transport  is  particularly  intense  in  the  swash  zone  and 
in  the  area  of  the  break  point  of  the  incident  waves.  However,  sand  is  in 
continual  motion  all  along  the  subaqueous  profile. 

13.  The  subaerial  portion  of  the  beach  also  changes.  Berms  build 
during  calmer  wave  conditions.  Sand  from  the  berm  may  be  blown  onshore  to 
form  dunes,  or  they  may  be  created  artificially.  If  the  water  level  rises, 
the  formerly  subaerial  portion  of  the  beach  will  be  unstable  relative  to 
sediment  moving  hydraulic  forces.  The  berm  and  dune  can  erode  by  overwash 
(transporting  sand  landward)  or  by  scour  and  collapse  (moving  beach  and  dune 
sediments  seaward). 

14.  Beach  profiles  are  typically  classified  into  two  basic  types, 
either  depositional  (also  normal,  summer,  swell,  step  type,  or  reflective)  or 
erosional  (also  storm,  winter,  bar  type,  or  dissipative)  as  shown  in  Figure 
3.  Other,  more  detailed  classifications  have  been  proposed  (Sunamura  and 
Horikawa  1975,  Wright  et  al.  1979,  Hattori  and  Kawamata  1981,  and  Wright  and 
Short  1983). 


Factors  Contributing  to  Profile  Change 

15.  Bruun  (1962),  Edelman  (1968,  1972),  Dean  (1977),  Vellinga  (1983b) 
and  others  have  identified  an  increase  in  water  level  as  the  most  important 
variable  relating  to  dune  and  shoreline  retreat.  During  a  storm,  the  "storm 
surge"  develops  as  a  result  of  the  superposition  of  the  astronomical  tide, 
wave  setup,  and  meteorological  (wind  and  pressure  induced)  surge.  This 
abnormally  high  water  level  permits  erosive  waves  to  attack  the  subaerial 
beach.  On  a  less  energetic  scale,  it  has  been  well  documented  (Inman  and 
Filloux  1960,  Otvos  1965,  and  Hattori  1983)  that  cross-shore  movement  of 
sediment  is  closely  related  to  the  tide  level. 

16.  A  change  in  water  level  does  not  cause  erosion;  wave  action  is 
required  to  move  the  sediment.  It  was  recognized  early  in  the  study  of 
coastal  processes  (Johnson  1952)  that  deepwater  wave  height  and  wave 
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Oppositional  Profile 


Figure  3.  Erosional  and  depositional  profile  shapes 
(after  Komar  1983) 

steepness  are  the  most  important  wave  parameters  controlling  beach  profile 
change.  Wave  steepness  is  defined  as  H0/L0  ,  in  which  H0  is  the  deepwater 
wave  height  and  L0  is  the  deepwater  wave  length.  In  linear  wave  theory, 

L0  =  gT^/(2^)  ,  in  which  g  is  the  acceleration  of  gravity,  and  T  is  the 
wave  period. 

17.  Johnson  found  that  the  critical  wave  steepness  for  delineating 
depositional  and  erosional  profiles  was  in  the  range  of  0.025  to  0.03. 

Higher  steepness  values  result  in  erosion;  lower  values  cause  deposition. 

For  example,  a  2-m-high,  8-sec-period  wave  (H0/L0  -  0.02)  implies  a 
depositional  profile;  whereas  a  2-m-high,  6-sec-period  wave  (H0/L0  -  0.036) 
results  in  an  erosional  profile.  In  this  example  the  wave  height  is 
constant,  and  the  period  is  different.  Combinations  of  wave  height  and 
period  are  not  completely  random.  Although  higher  waves  usually  have  longer 
periods,  storm  waves  typically  change  from  high  steepness,  short-period  waves 
to  lower  steepness,  long-period  waves  as  the  storm  passes. 

18.  Other  factors  affecting  profile  change  include  the  beach  grain 
size  (Iwagaki  and  Noda  1963)  or  fall  velocity  (Dean  1973),  mean  beach  slope 
(e.g.,  Sunamura  and  Horikawa  1975,  and  Hattori  and  Kawamata  1981),  and 
magnitude  of  the  wave  height  (Saville  1957).  Hallermeier  (1984)  has  recently 
presented  a  classification  of  the  occurrence  of  nearshore  deposition  and 
erosion  through  use  of  a  sediment  mobility  number  and  the  Ursell  number. 

Each  of  the  above  cited  papers  has  examined  or  presented  one  or  more  criteria 
for  predicting  the  tendency  to  form  either  an  erosional  or  depositional 
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Erosion 


19.  Starting  from  a  depositional  profile  characterized  by  a  high,  wide 
berm,  steep  foreshore  slope,  and  a  relatively  flat  offshore,  if  the  incident 
wave  steepness  exceeds  the  critical  steepness,  the  beach  begins  to  erode.  If 
the  surge  is  high  enough,  storm  waves  move  sand  from  the  berm  and  dune  off¬ 
shore  in  the  region  shoreward  of  the  main  wave  breaking  point.  As  the 
material  moves  offshore,  the  beach  slope  in  the  surf  zone  becomes  milder.  As 
a  result,  one  or  more  offshore  bars  may  form,  and  existing  bars  may  move 
farther  seaward  (Birkemeier  1985a).  Figure  4  illustrates  the  offshore  move¬ 
ment  of  a  bar  resulting  from  the  passage  of  a  series  of  three  storms  at  Duck, 
North  Carolina.  The  last  and  largest  storm,  occurring  13-15  November  1981, 
caused  the  most  beach  erosion  and  greatest  seaward  shift  of  material. 


Figure  4.  Offshore  bar  movement  resulting  from 
three  storms  at  Duck,  North  Carolina 
(Birkemeier  1985a) 


20.  The  mechanisms  controlling  offshore  bar  formation  and  movement  are 
not  well  understood.  Two  main  mechanisms  are  known  to  exist.  The  classical 
mechanism  is  that  the  bars  are  formed  at  the  break  point  of  storm  waves;  the 
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other  mechanism  is  that  bars  are  formed  by  the  convergences  and  divergences 
of  mean  flows  associated  with  partially  standing  waves  with  long  periods 
(Bowen  and  Inman  1971,  Short  1975,  and  Holman  and  Bowen  1982).  What  is 
understood  is  that  bars  provide  a  natural  protection  to  the  beach  by  causing 
incoming  waves  to  break  farther  offshore,  reducing  the  energy  reaching  the 
beach.  Offshore  bars  also  provide  a  mechanism  for  the  temporary  storage  of 
eroded  sediments  as  part  of  the  cycle  toward  beach  recovery. 

Beach  recovery 

21.  As  a  storm  subsides,  the  wave  height  decreases,  and  the  wave 
period  usually  increases,  eventually  resulting  in  a  wave  steepness  below  the 
critical  value  for  erosion.  Then  the  beach  recovery  process  begins.  In  the 
inner  part  of  the  surf  zone,  sand  is  transported  onshore,  and  a  new  berm  is 
created  by  the  runup.  Because  sand  is  removed  from  the  surf  zone  and  carried 
onshore,  the  slope  of  the  surf  zone  becomes  steeper.  The  process  of  berm 
creation  has  been  little  studied  and  is,  therefore,  not  well  understood.  The 
rate  of  beach  recovery  can  be  quite  rapid;  for  example,  Birkemeier  (1979) 
reported  that  within  one  day  of  the  peak  of  a  storm,  over  51  percent  of  the 
eroded  material  had  returned  to  the  beach  at  Long  Beach  Island,  New  Jersey. 

22.  In  the  offshore,  the  storm- induced  bar  will  move  landward  and 
ultimately  onto  the  beach.  The  rate  of  recovery  can  vary  considerably 
depending  on  wave  conditions  and  the  depth  over  the  bar.  Figure  5 


Figure  5.  Slow  recovery  of  an  offshore  bar  at  Duck 
North  Carolina  (Birkemeier  1985a) 


illustrates  the  slow  recovery  of  the  bar  shown  in  Figure  4.  Recovery  did  not 
begin  until  the  following  February  and  continued  for  the  next  6  months,  an 
extended  period  of  low  wave  height  conditions.  Recovery  of  shallower  bars 
can  be  rapid.  Sallenger,  Holman,  and  Birkemeier  (1985)  measured  a  rate  of 
onshore  bar  migration  within  the  surf  zone  of  1.2  m/hr  during  the  waning 
stages  of  a  storm  when  the  waves  were  2  m  high  with  a  14-sec  period. 

23.  Given  a  sufficiently  long  period  (several  months  or  more)  of  low 
waves,  the  storm  bar  may  disappear  (e.g.,  Birkemeier  1985a).  If,  however, 
the  storm  bar  is  located  in  very  deep  water,  the  bar  may  nor  be  erased  over 
an  extended  time  scale  of  many  years.  In  this  case,  if  no  other  source  of 
sand  is  available,  the  beach  suffers  from  an  apparent  permanent  erosion  (Dean 
1976)  . 


Overwash 


24.  "Overwash"  is  a  process  which  occurs  when  the  storm  surge  and  wave 
uprush  overtop  the  foredune  line  causing  water  and  sand  to  move  inland. 
Overwash  is  of  major  importance  during  extreme  storms  when  a  significant 
percentage  of  a  duneline  may  be  overtopped  or  breached  (Figure  6).  It  is  an 
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Figure  6.  Breached  and  overtopped  dune  on  Grand  Isle,  Louisiana, 
following  Hurricane  Juan  in  1985 
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important  geological  process  since  it  provides  a  mechanism  for  barrier 
islands  to  move  and  to  increase  in  elevation.  Correct  identification  of 
potential  overwash  areas  is  also  important  in  the  determination  of  V-zones. 

25.  Most  overwash  research  has  been  geologic  in  nature,  with  few 
engineering  studies  having  been  made.  A  comprehensive  evaluation  of  the 
literature  was  conducted  and,  to  the  authors'  knowledge,  no  analytical  model 
of  the  overwash  process  exists  except  that  developed  in  an  M.  S.  thesis  by 
Williams  (1978).  However,  the  laboratory  tests  and  theory  of  Williams  were 
inconclusive.  Neither  the  two-dimensional  (2-D)  dune  response  to  overtopping 
nor  the  three-dimensional  problem  associated  with  dune  breaching  has  been 
adequately  studied  to  permit  numerical  simulations  of  the  overwash  process. 
The  Multiple  Shore-Breaking  Wave  Transformation  (MSBWT)  model  of  Balsillie 
(1984c,  1985b)  has  provision  for  determining  wave  heights  over  both  breached 
and  flooded  profiles  as  well  as  a  mechanism  for  estimating  overwash 
deposition,  but  the  model  is  not  fully  operational  or  verified.*  A  detailed 
discussion  on  this  model  is  given  in  Part  III.  Other  studies  quantify  the 
percentage  of  island  overwashed  (Pierce  1969,  and  Boc  and  Langfelder  1977)  or 
address  volumetric  changes  (Schwartz  1975,  and  Leatherman  1976,  1981).  The 
results  of  the  literature  survey,  together  with  a  summary  of  overwash 
knowledge  for  specific  case  examples,  are  contained  in  Appendix  A.  A 
procedure  for  accounting  for  overwash  in  V-zone  determination  is  discussed  in 
Part  VI. 


Equilibrium  Beach  Profile  Concept 


26.  It  is  evident  to  the  most  casual  observer  that  the  beach  profile 
has  an  apparent  constancy  of  shape.  When  we  visit  a  particular  beach,  the 
profile  looks  more  or  less  the  same,  year  after  year.  This  is  most  easily 
explained  using  the  concept  of  an  "equilibrium  profile"  which  continually 
evolves  to  be  in  equilibrium  with  whatever  wave  and  water  level  conditions 
exist.  This  concept  is  fundamental  to  existing  predictive  models  of  dune 
erosion. 


Information  obtained  from  personal  communication  with  J.  H.  Balsillie, 
Division  of  Beaches  and  Shores,  Florida  Department  of  Natural  Resources, 
Tallahassee,  Fla.,  June  1986. 


27.  On  an  idealized  2-D  profile,  waves  provide  the  mechanism  for 
cross-shore  movement  of  sand,  with  movement  either  onshore  or  offshore 
depending  on  wave  steepness.  If  constant  wave  and  water  level  conditions 
continue  for  an  extended  period  of  time,  an  equilibrium  shape  will  be 
approached.  The  development  of  this  equilibrium  profile  has  been  demon¬ 
strated  with  physical  model  tests  using  monochromatic  waves  (e.g.,  Saville 
1957,  Swart  1974,  and  Chesnutt  1978).  In  general,  laboratory  experiments 
have  used  constant  water  depths  and  an  initial  plane  beach  profile.  Measured 
changes  are  greatest  at  the  start  of  the  test,  gradually  decreasing  with 
time.  Since  true  equilibrium  is  difficult  to  obtain  in  the  laboratory,  it  is 
generally  believed  to  be  seldom  reached  in  nature. 

28.  Profile  changes  resulting  from  changes  in  wave  characteristics  may 
be  thought  of  as  perturbations  on  a  smoother,  more  idealized  profile  shape 
that  more  or  less  retains  a  constant  form  in  equilibrium  with  the  mean  water 
level,  mean  waves,  and  beach  grain  size.  As  the  water  level  changes,  this 
equilibrium  profile  shifts  accordingly,  moving  inland  and  upward  with 
increasing  water  level.  The  shape  of  the  profile  is  maintained  by  erosion  of 
the  beach  and  deposition  offshore.  Profiles  are  more  likely  to  be  in 
equilibrium  with  long-term  effects  such  as  sea  level  rise  but  not  with  short- 
duration  storm  surges. 

29.  In  an  empirical  study  of  beaches  on  the  north  Danish  coast  and  in 
California,  Bruun  (1954)  found  that  the  representative,  gradually  deepening 
beach  profile  was  well  described  by  a  simple  equation  of  the  form 

h(x)  =  Ax^/3  (1) 

where 

h  =  depth 

x  =  distance  from  the  shoreline 

A  =  dimensional  factor  having  units  of  length  to  the  one  third  power 

A  definition  sketch  is  given  in  Figure  7.  Dean  (1977)  found  the  same 
relationship  empirically  for  502  profiles  along  the  Atlantic  Ocean  coast  of 
the  United  States.  Subsequent  work  (e.g.,  Hughes  1978,  and  Moore  1982)  has 
verified  the  basic  validity  of  Equation  1  for  beach  profiles  in  various  parts 
of  the  world  and  for  profiles  generated  in  laboratory  flumes  and  basins  with 
sand  bottoms.  Although  the  two-third  power  law  does  not  exhibit  bars  and 
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troughs,  it  represents  a  best  fit  description  of  a  profile  passing  through 
such  features. 


Figure  7.  Definition  sketch  for  the  equilibrium  beach  profile 

(after  Kriebel  1982) 

30.  Dean  (1977)  gave  a  plausible  derivation  leading  to  Equation  1.  He 
showed  that  if  energy  dissipation  per  unit  volume  of  water  in  the  surf  zone 
were  assumed  to  be  the  dominant  factor  controlling  profile  shape,  the  two- 
third  power  law  would  result.  In  addition,  the  general  functional  dependence 
of  the  shape  factor  A  on  the  dissipation  could  be  predicted.  Dean  (1977), 
Hughes  (1978),  and  Moore  (1982)  have  demonstrated  that  A  depends  in  a 
rational  way  on  the  grain  size  or  fall  velocity  of  the  beach  material.  A 
grain  size  of  0.25  mm  corresponds  to  a  value  of  A  of  approximately  0.13 
ml/3  (Moore  1982).  For  larger  grain  sizes,  A  increases  producing  a  steeper 
profile,  as  observed  in  nature.  Figure  8  shows  the  dependency  of  A  on  the 
grain  size  according  to  Moore  (1982). 

31.  In  the  derivation  of  Dean  (1977),  the  required  energy  dissipation 
is  calculated  from  an  assumed  linear  and  constant  wave  height  decay  with 
depth,  using  small  amplitude  wave  theory.  This  procedure  automatically 
restricts  applicability  of  the  method  to  a  surf  zone  with  spilling  breakers 
for  which  the  breaker  height  and  water  depth  are  linearly  related.  However, 
if  a  smoothing  mechanism  is  posited  which  slowly  shifts  sand  along  the 
profile,  it  can  be  assumed  that  the  equilibrium  profile  will  extend  to  at 
least  the  depth  corresponding  to  the  depth  of  the  greatest  breaking  waves 
during  storms.  The  smoothing  mechanism  might  simply  be  the  back-and-f orth 
sand  movement  associated  with  the  wave  orbital  velocities. 
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Sediment  Diameter  (mm) 

Figure  8.  Shape  factor  A  versus  sediment  diameter 

(Moore  1982) 

32.  The  equilibrium  profile  concept  is  useful  for  two  purposes:  it 
provides  a  convenient  means  of  visualizing  the  process  of  profile  change; 
and,  through  Equation  1,  it  provides  a  reasonable  description  of  the  general 
shape  of  the  beach  profile. 

Natural  Variation  in  Measured  Profiles 

33.  As  already  mentioned,  the  beach  profile  is  in  constant  motion 
under  the  imposed  hydraulic  forces  which  move  sediment  across  and  along  the 
shore.  The  profile  thereby  undergoes  notable  short-term  macroscale  changes 
(order  of  days  and  tens  of  cubic  meters)  apart  from  seasonal  and  long-term 
changes.  On  real  and  laboratory  beaches,  these  changes  do  not  have  perfect 
symmetry  alongshore;  there  is  always  some  degree  of  three-dimensionality. 
This  longshore  variability  may  result  from  the  presence  of  shore  protection 
structures  such  as  jetties  or  groins,  or  the  variation  may  be  naturally 
generated  by  the  hydraulic  conditions  (e.g.,  by  beach  cusps  or  rip  currents) 
Varying  geological  characteristics  along  the  coast  will  also  produce  markedl 
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different  longshore  conditions.  In  any  case,  we  must  quantify  the  naturally 
occurring  variability  in  order  to  establish  accuracy  requirements  for  a  dune 
erosion  mode  1 . 

34.  Birkemeier,  Savage,  and  Leffler  (in  preparation)  present  a  large 
collection  of  storm  erosion  data,  46  cases  from  seven  different  beach 
localities  and  13  different  storm  events  (Table  1).  Included  in  the  study  is 
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Locality 

Tide  gage  used 

Note:  a  M-"  Indicates  irtcnown  <tota. 

*  -  peak  sirge,  maasired  above  M3W). 

NAU5T 

touset  Beacti,  Cape  Cod,  toss. 

Boston  tortcr ,  Mass. 

M  -  hlndcasted  wave  hel^rts  at  9-m  depth. 

MI  SO 

Mlstjjanlcut  Beacti,  R.I. 

tontatJc  Pt.,  N.V. 

f  -  sirga  data  from  Sandy  took,  N.J. 

VCST 

WBSthanpton  Beacti,  N.V. 

tontaiJc  Pt.,  N.V. 

ft  -  based  on  visual  breaking  wave  heltfrt. 

JOES 

Jones  Beach,  N.V. 

Sandy  took,  N.J. 

$  -  gage  masirenent  In  10-m  water  depth. 

LBI 

Long  Beach  Island,  N.J. 

Atlantic  City,  N.J. 

ATLANT 

Atlantic  City,  N.J. 

Atlantic  City,  N.J. 

ILDLAM 

ludlaw  Beach,  N.J. 

Atlantic  City,  N.J. 

a  total  of  549  different  measured  profile  changes.  Though  no  hurricanes  or 
100-year  events  are  included,  it  is  a  unique  data  set  with  which  to  examine 
natural  variability  in  profile  change. 
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35.  Each  case  was  specifically  selected  from  a  larger  group  of  surveys 
to  isolate  the  effects  of  single  storms.  Because  all  poststorm  surveys  were 
conducted  within  9  days  of  the  storm,  the  data  contain  minimal  effects  of 
poststorm  beach  recovery.  For  each  case,  shoreline  and  volume  changes 
relative  to  mean  sea  level  (MSL)  are  computed  and  presented  for  a  number  of 

JL 

beach  profile  lines.  Birkemeier,  Savage,  and  Leffler  (in  preparation)  also 
report  volume  changes  based  on  half-meter  contours  of  elevation.  Because  the 
present  study  is  concerned  with  erosion  above  the  surge  level  (including  the 
astronomical  tide),  these  data  were  used  to  compute  the  volume  changes  above 
the  peak  measured  water  level. 

36.  Also  included  in  Table  1  is  the  wave  height  and  water  level  infor¬ 
mation  associated  with  each  storm.  The  wave  data  were  hindcast  according  to 
Jensen  (1983).  Water  level  measurements  were  obtained  from  nearby  tide 
gages.  The  peak  levels  given  in  Table  1  were  computed  relative  to  MSL.  Note 
that  because  the  water  level  data  were  not  obtained  at  each  beach,  a  true 
measure  of  the  local  surge  and  the  volume  change  above  the  surge  level  is  not 
available. 

Shoreline  position 

37.  A  usual  measurement  of  profile  change  is  the  change  in  shoreline 
position,  which  is  defined  as  the  MSL  intercept.  Birkemeier,  Savage,  and 
Leffler  (in  preparation)  have  found  that  generally  both  the  MSL  shoreline 
position  and  the  beach  slope  at  MSL  are  relatively  insensitive  to  storm 
changes,  regardless  of  locality  or  storm.  For  example,  of  the  549  profile- 
survey  combinations  for  which  an  MSL  shoreline  change  could  be  computed,  81 
percent  eroded  based  on  volumetric  changes  (regardless  of  magnitude);  whereas 
only  54  percent  of  the  lines  had  erosional  shorelines.  Although  poststorm 
recovery  (which  would  affect  the  shoreline  position)  is  certainly  a  factor  on 
the  46  percent  of  the  lines  which  showed  shoreline  accretion,  there  is 


considerable  evidence  indicating  that  the  MSL  shoreline  position  effectively 
acts  as  a  pivot  point. 

38.  Figure  9  shows  the  distribution  of  median  shoreline  changes  by 
storm  and  locality.*  For  each  storm  and  locality,  a  small  "box  plot"  illus¬ 
trates  the  distribution  of  measured  profile  changes.  The  portion  of  the  box 
above  or  below  the  dashed  line,  respectively,  indicates  either  accretion  or 
erosion.  Overall,  the  average  median  shoreline  change  was  small,  only  -0.9 
m,  with  a  mean  hinge  range  (difference  of  the  hinges)  of  4.8  m  (+1.7  to  -3.1 
m) .  Median  shoreline  position  changes  with  an  absolute  value  less  than  2  m 
were  recorded  for  39  percent  of  the  cases.  Only  8  percent  of  the  shoreline 
changes  had  an  absolute  value  greater  than  10  m.  It  is  interesting  to  note 
that  the  range  of  variation  between  cases  (indicated  by  the  heights  of  the 
boxes  in  Figure  9)  is  relatively  small  and,  with  a  few  exceptions,  measured 
variations  are  as  similar  between  storms  at  one  locality  as  they  are  between 
localities.  Although  shoreline  position  is  traditionally  used  for  computing 
long-term  erosion  rates,  based  on  this  data  set  use  of  the  change  in 
shoreline  position  for  quantifying  storm  erosion  appears  to  be  limited. 

Volume  changes 

39.  Figure  10  plots  the  distribution  of  volume  change  above  MSL  for 
each  locality  and  storm.  Unlike  the  shoreline  changes,  the  volume  changes 
show  more  consistent  erosion  and  more  variation  between  storms.  Misquamicut 
and  Ludlam  Beach  had  the  smallest  ranges  in  variation  between  profile  lines; 
whereas  Nauset  Beach  and  Jones  Beach  had  relatively  large  ranges.  Individual 
profile  change,  represented  by  the  extreme  values,  can  be  quite  large,  up  to 
-150  m^/m  for  one  Atlantic  City,  New  Jersey  profile.  Though  it  is  difficult 
to  intercompare  storms,  the  data  in  Figure  10  illustrate  the  large  amount  of 
variation  which  naturally  occurs  both  along  a  beach  and  between  beaches  sub¬ 
jected  to  the  same  storm. 

40.  In  order  to  account  for  this  natural  variation  caused  by  storms, 
the  data  shown  in  Figure  10  were  used  to  comp'ute  a  "variability  factor"  or 

*  The  median  value  is  preferred  over  the  mean  of  a  number  of  profile  lines 
because  it  is  insensitive  to  single  profiles  with  extreme  changes  and  should 
be  more  representative  of  the  overall  locality  change.  Plotted  changes  show 
both  upper  and  lower  extreme  values  and  "hinge"  values.  Hinge  values  are 
defined  as  the  25  and  75  percentiles  for  the  profiles  at  each  locality. 


multiplier  to  be  used  with  a  dune  erosion  estimate.  Assuming  that  the 
estimate  was  close  to  the  median  volume  change  on  a  beach  and  that  a 
desirable  prediction  should  cover  75  percent  of  the  changes  (or  any 
acceptable  percentage),  then  an  average  multiplier  of  the  median  value  could 
be  computed  to  predict  the  75  percent  change.  Of  the  45  cases  reported  by 
Birkemeier,  Savage,  and  Leffler  (in  preparation)  which  had  median  above-MSL 
volume  losses  in  excess  of  -1  m^/m  of  erosion,  a  median  variability  factor  of 
1.6  was  computed  and  is  plotted  along  with  the  distribution  of  factors  in 
Figure  11.  Also  shown  in  Figure  11  is  a  similar  distribution  of  factors 
based  on  the  computed  volume  changes  above  the  surge  level.  Only  35  of  the 
46  cases  given  in  Table  1  had  median  above-surge  level  volume  losses  in 
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Figure  11.  Distribution  of  computed  variability  factors 

excess  of  -1  m^/m.  It  is  interesting  to  note  that  the  median  above-surge 
variability  factor  of  1.5  is  nearly  identical  to  the  above-MSL  factor.  Use 
of  the  75  percent  factor  of  2.1  is  recommended.  This  value  can  be  inter¬ 
preted  as  the  multiplier  of  the  median  required  to  include  volumetric  erosion 
expected  on  75  percent  of  the  shoreline.  A  higher  value,  for  example  2.5, 
would  increase  the  percentage  included  and  would  conservatively  account  for 
the  uncertainty  in  the  dune  erosion  estimate.  This  specific  value  of  the 
variability  factor  should  be  viewed  as  preliminary,  to  be  revised  as  more 
data  become  available. 


41.  In  principle,  the  value  of  the  variability  factor  could  be 
dependent  on  surge  height,  wave  height,  and  the  median  change.  No  such 
dependency  was  identified  in  the  data  analyzed. 

42.  The  need  for  a  variability  factor  has  significant  implications 
relative  to  the  required  accuracy  of  the  method.  Obviously,  with  such  a  high 
degree  of  variability,  a  precise  estimate  of  erosion  along  an  individual 
profile  is  not  required;  however,  a  reasonable  estimate  of  the  maximum  (here 
chosen  as  the  75th  percentile)  change  which  may  occur  on  any  profile  is 
required.  An  important  added  consideration  is  the  large  degree  of  uncer¬ 
tainty  concerning  the  long-term  stability  of  a  particular  beach  prior  to  the 
arrival  of  the  100-year  storm. 

43.  Although  this  report  does  not  address  additional  requirements 
introduced  by  long-term  erosion  (over  periods  of  many  years),  the  process 
cannot  be  ignored.  Birkemeier  (1979)  reported  on  significant  changes 
resulting  from  only  a  2-year  return  period  storm.  The  profile  shape  at  any 
particular  time  is  related  to  the  integrated  effect  of  all  previous  storm 
(and  nonstorra)  waves.  Since  predictions  are  invariably  based  on  a  limited 
amount  of  beach  survey  data,  the  accuracy  and  usefulness  of  the  prediction 
will  decrease  with  time.  As  an  extreme  example,  a  dune  sufficiently  sub¬ 
stantial  to  survive  a  major  storm  may  disappear  under  the  combined  attack  of 
a  number  of  small  storms  over  a  number  of  years.  The  Grand  Isle  case 
discussed  in  Part  I  is  just  one  example  of  a  scenario  repeated  along  almost 
every  shoreline  which  is  undergoing  long-term  erosion. 


PART  Ills  PREDICTIVE  MODELS  OF  PROFILE  CHANGE  AND  DUNE  EROSION 


44.  This  section  identifies  the  criteria  by  which  the  applicability 
and  suitability  of  a  numerical  model  for  simulating  profile  change  can  be 
judged.  The  criteria  define  the  properties  of  an  ideal  model,  one  that 
neither  presently  exists  nor  is  expected  to  be  developed  in  the  immediate 
future.  The  list  is  given  to  provide  a  standard  or  checklist  with  which 
existing  profile  change  models  can  be  evaluated. 

45.  A  dune  erosion  model  need  not  possess  all  the  properties  of  a 
profile  change  model  because  a  dune  erosion  model  is  aimed  at  simulating  an 
extreme  erosion  event.  The  process  of  beach  recovery  is  assumed  to  be 
secondary  and  can  be  neglected  at  the  first  stage  of  dune  erosion  model 
development.  Moreover,  as  discussed  earlier,  the  accuracy  requirements  of  a 
model  for  V-zone  determination  are  much  less  restrictive  than  those  required 
for  an  empirical  process-response  numerical  sediment  transport  model. 
Therefore,  a  subset  of  the  criteria  defining  the  properties  of  an  ideal  beach 
profile  model  is  identified  for  use  in  evaluating  dune  erosion  numerical 
models . 

46.  Several  of  the  more  recent  theoretical  models  of  cross-shore 
sediment  transport  and  beach  profile  change  are  briefly  described.  Then  a 
more  detailed  evaluation  of  one  of  the  models  is  made  from  the  perspective  of 
potential  use  as  a  dune  erosion  model. 

Required  Characteristics  of  Beach  Profile  Change  Models 

47.  An  attempt  was  made  to  compile  a  complete  list  of  properties  and 
characteristics  required  for  an  ideal  numerical  model  of  profile  change,  but 
it  is  recognized  that  much  has  yet  to  be  learned  about  the  processes  of  wave 
deformation  in  the  surf  zone,  beach  profile  change,  dune  erosion,  and 
overwash.  The  individual  items  or  requirements  listed  below  are  not 
completely  independent  of  each  other  but  are  presented  separately  for 
clarity.  These  14  items  have  been  placed  into  four  groups:  input  data  and 
initial  and  boundary  conditions,  calculation  procedures,  model  properties, 
and  extendabi 1 ity .  Other  combinations  of  properties  are  possible  also.  For 
example,  Dally  (1980)  and  Dally  and  Dean  (1984)  have  presented  five  criteria 
which  a  "good"  beach  profile  model  should  satisfy.  The  four  criteria  for  an 


ideal  beach  profile  change  model  discussed  in  this  report  include  the 
following : 


f 


a .  Input  data,  initial  conditions,  and  boundary  conditions. 

(1)  Deepwater  wave  height  and  wave  period  (entering  in  part 
as  the  wave  steepness)  should  be  among  the  basic  input 
variables  to  drive  the  model  or  to  establish  an  offshore 
wave  boundary  condition. 

(2)  The  properties  of  the  beach  material  (fall  velocity  or 
representative  grain  size  and  composition)  must  be 
included.  It  may  be  necessary  to  allow  for  natural 
sorting  along  the  profile  in  some  cases. 

(3)  The  profile  shape  is  required  as  an  initial  condition. 

(4)  An  offshore  boundary  condition,  specifying  the  maximum 
depth  to  which  sand  is  transported  across  shore  in 
significant  amounts,  is  required. 

(5)  An  onshore  boundary  condition  is  required.  This  boundary 
condition  must  allow  for  berm  and  dune  erosion,  as  well 
as  berm  buildup. 


b. 


Calculation  procedures. 

(1) 


Calculation  procedures  for  the  local  incident  waves, 
infragravity  waves,  and  wave-induced  currents  (orbital 
velocities,  mass  transport,  or  undertow)  in  the  region 
between  the  offshore  and  onshore  boundary  points  are 
required.  This  calculation  includes  breaking  wave 
criterion  and  broken  (surf  zone)  wave  model. 


(2) 


The  model  must  allow  for  calculation  or  specification  of 
the  mean  water  level,  including  the  effects  of  astronomi¬ 
cal  tide,  meteorologi cal  surge,  wave  setup,  and  runup. 


(3) 


Calculation  procedures  for  cross-shore  sediment  transport 
rates  are  necessary.  These  procedures  would  include 
threshold  effects  and  encompass  spatially  and  temporally 


(4) 
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c^.  Model  properties . 


(1)  Formation,  movement,  and  erasure  of  bars  must  be 
represented . 
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(2)  If  the  wave  conditions  and  water  level  are  held  constant, 
the  model  should  allow  spatially  and  temporally  accurate 
calculation  of  both  erosional  and  depositional  profiles 
and  produce  an  equilibrium  profile  shape. 

(3)  The  model  must  be  verified  to  produce  accurate  results 
for  the  range  of  conditions  for  which  it  is  expected  to 
be  applied. 

d .  Extendability . 

(1)  Since  a  deterministic  model  is  not  expected  to  be 
practicable  for  field  use,  the  developed  model  must  be 
compatible  with  stochastic  or  probabilistic  prediction 
procedures . 

(2)  Although  as  a  first  step  a  monochromatic  wave  input  is 
expected  to  be  adequate,  refinements  should  include  the 
randomness  of  wave  height,  direction,  and  period. 


48.  The  deepwater  wave  height  and  wave  period  are  usually  available  or 
may  be  estimated  with  acceptable  accuracy  (item  a(l)).  It  is  difficult  to 
specify  the  beach  material  grain  size  along  the  profile  as  required  by  a(2). 
The  minimum  requirement  is  knowledge  of  an  average  grain  size  or  fall 
velocity  for  the  target  beach.  For  a(3),  in  principle,  the  exact  initial 
profile  is  required.  This  information  will  not  be  available  in  a  practical 
situation.  Either  a  historic  profile  or  a  representative  profile  shape  will 
have  to  be  used  in  applications. 

49.  For  the  offshore  sediment  transport  boundary  condition  (item 
a(4)),  the  work  of  Hallermeier  (1979),  as  modified  by  Birkemeier  (1985b),  is 
probably  acceptable  at  the  first  stage  of  modeling.  At  present,  little 
knowledge  exists  to  formulate  a  first-principles  onshore  boundary  condition 
(item  a(5)).  A  rigorous  onshore  boundary  condition  will  probably  incorporate 
a  number  of  boundary  conditions  for  the  foreshore,  for  the  berm,  and  for  the 
dune.  For  the  present,  simple  geometrical  arguments  may  have  to  be  used. 

50.  For  item  b(l),  at  present,  accurate  calculation  of  the 
cross-shore  current  flow  through  the  water  column  in  the  surf  zone  and 
seaward  of  the  breakers  is  not  possible.  Water  levels  (item  b(2))  can  be 
estimated  with  sufficient  accuracy  for  the  purpose  of  a  dune  erosion  model. 
The  requirement  of  valid  cross-shore  sediment  transport  rate  relationships 
for  all  regions  along  the  profile  (item  b(3))  cannot  be  satisfied  by  any 
single  formulation  or  combination  thereof  presently  available.  Beach  profile 
change  and  the  sediment  transport  rate  are  connected  through  a  continuity 


equation  for  beach  material,  and  the  equation  must  be  properly  expressed  in 
numerical  form  (item  b(4)). 

51.  Requirement  c(l)  states  that  the  numerical  formulation  must  be 
capable  of  representing  all  important  morphological  features.  The  scheme  must 
produce  a  profile  which  is  in  equilibrium  with  imposed  waves  of  constant 
properties  (c(2)).  This  requirement  is  important  because  it  is  often 
possible  for  a  numerical  scheme  to  appear  to  be  stable  for  a  small  number  of 
calculation  time-steps  but  will  diverge  if  allowed  to  continue  for  a  longer 
simulation  interval.  This  situation  implies  that  the  intermediate  computed 
results  may  be  spurious.  Naturally,  any  model  must  be  calibrated  and 
verified  by  suitable  field  data  (item  c(3)).  Only  then  can  it  be  confidently 
used  for  predictions. 

52.  Concerning  the  requirement  for  extension  to  a  probabilistic 
description  (item  d(l)),  it  has  been  shown  in  Part  II  that  a  relatively  wide 
range  in  natural  variability  exists  in  beach  profile  change.  In  addition,  an 
artificial  variability  exists  because  the  initial  conditions  required  to  set 
up  a  numerical  model  of  profile  change  (wave  conditions,  initial  profile 
shape,  etc.)  are  not  accurately  known.  The  most  logical  approach  would  be  to 
calculate  a  range  of  values  for  profile  change,  within  which  the  actual 
change  would  be  expected  to  lie.  This  approach  leads  to  a  stochastic  or 
probabilistic  description.  For  item  d(2),  it  is  known  that  random  waves 
exhibit  somewhat  different  properties  from  those  of  monochromatic  waves.  At 

a  later  stage  in  modeling,  the  random  nature  of  the  wave  field  should  be 
taken  into  account  as  well  as  sediment  transport  under  random  wave  action. 


Required  Characteristics  of  Dune  Erosion  Models 


53.  A  pragmatic  engineering  simulation  model  for  dune  erosion  need  not 
possess  all  the  characteristics  of  a  theoretically  ideal  beach  profile  change 
model.  In  this  subsection,  a  subset  of  parameters,  expected  to  comprise  a 
first-stage  dune  erosion  model,  is  listed. 

a^.  Offshore  wave  height  and  wave  period.  (Local  wave  height  may 
or  may  not  be  necessary.) 

b.  Water  level,  including  astronomical  tide,  meteorological  tide, 
and  wave  setup . 


Representative  grain  size  or  fall  velocity  for  the  profile. 

d.  Representative  initial  shape  of  the  profile,  including  the 
dune . 

e.  Onshore  boundary  condition(s)  for  dune  erosion. 

f_.  Phenomenological  calculation  procedure  for  the  cross-shore 

sediment  transport  rate  or,  alternatively,  direct  calculation 
of  the  morphologic  change  of  the  profile  based  on  a  given 
idealized  form. 

Calculation  procedure  for  computing  dune  erosion  that  contains 
a  temporal  dependence. 

h.  Verification  with  field  data  over  the  range  of  conditions  for 
which  the  model  will  be  applied. 


The  conditions  for  defining  a  pragmatic  dune  erosion  model  are  fewer  and 
considerably  less  rigorous  compared  to  those  defining  an  ideal  beach  profile 
change  model. 


Numerical  Models  of  Cross-Shore  Sediment  Transport 
and  Beach  Profile  Change 


54.  This  subsection  introduces  several  recent  models  which  might  be 
considered  as  candidates  for  developing  a  "first  principles"  model  of  dune 
erosion.  It  will  be  concluded  that  none  of  the  models  can  be  readily  adapted 
for  immediate  use  as  a  dune  erosion  model. 

55.  The  models  may  be  classified  into  one  of  two  groups  according  to 
whether  they  rely  on  basic  sediment  transport  equations  or  on  empirical 
results.  Models  developed  from  relations  for  the  cross-shore  sediment 
transport  rate  include  those  of  Bailard  (1981,  1982,  1983,  and  1985),  Bowen 
(1980),  Dally  and  Dean  (1984),  Leont'ev  (1985),  Moore  (1982),  and  Yang 
(1981).  Empirically  based  models  include  those  of  Swart  (1974,  1975,  and 
1977),  Kajima  et  al.  (1983),  and  Swain  and  Houston  (1983,  1985).  Three 
models  developed  specifically  for  estimating  dune  erosion  and  profile  change 
are  discussed  in  Part  IV. 

Theoretical  models 

56.  The  transport  rate  equations  of  Bowen  (1980)  and  Bailard  (1981, 
1982,  1983,  and  1985)  are  based  on  the  energetics  approach  of  Bagnold  (1963, 
1966).  Both  contain  provisions  to  calculate  bed  load  and  suspended  load 
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sediment  transport.  Bowen's  model  has  not  been  developed  further  than  the 
initial  formulation.  Bailard  (1985)  states  that  the  performance  of  his  model 
was  "disappointing"  in  a  comparison  of  predicted  and  measured  profile  change. 
Bailard's  work  is  described  in  the  next  subsection.  Dally  and  Dean  (1982) 
developed  a  beach  profile  numerical  model  by  assuming  mean  wave-induced  flows 
are  responsible  for  the  sediment  transport.  Only  suspended  load  was 
considered.  Leont'ev  (1985)  uses  an  energetics-based  approach  and  assumes 
bed-load  transport  is  directed  solely  onshore;  whereas  suspended  load  is 
directed  solely  offshore.  Mean  flows  are  calculated.  Yang's  (1981)  model 
employs  both  mean  and  fluctuating  currents,  and  the  suspended  load  (no  bed 
load  is  included)  is  assumed  to  be  transported  only  offshore.  Moore  (1982) 
uses  an  equilibrium  profile-based  transport  model  plus  a  surf  zone  wave 
model . 

57.  It  is  remarkable  that  despite  quite  different  assumptions,  all  of 
the  above  mentioned  models  were  demonstrated  by  the  authors  to  represent 
beach  profile  change  to  some  extent.  Notably,  in  the  development  of  the 
cited  models,  little  discussion  was  given  to  the  onshore  and  offshore 
boundary  conditions  and  to  the  approach  to  equilibrium.  Also,  only  very 
limited  verification  data  were  employed.  In  particular,  no  comparison  was 
made  to  data  on  the  sediment  transport  rate  but  only  to  profile  change.  In 
order  to  further  explore  these  models,  it  is  useful  to  examine  one  of  them, 
Bailard's  sediment  transport  model,  in  detail. 

Sediment  transport  model  of  Bailard 

58.  Bailard  (1981)  derived  equations  for  both  longshore  and  cross¬ 
shore  sediment  transport  on  an  assumed  idealized  beach  of  constant  slope. 

The  model  is  based  on  the  energetics  approach  of  Bagnold  (1963,  1966).  Bowen 

(1980)  presented  a  similar  model.  Since  the  model  of  Bailard  has  undergone 
refinement  and  testing  (Bailard  1982,  1983,  and  1985),  Bailard's  work  was 
selected  for  closer  examination,  and  discussion  will  be  limited  to  the  cross¬ 
shore  component.  It  should  be  mentioned  that  Bailard's  work  mainly  concerns 
the  theoretical  derivation  and  interpretation  of  a  general  sediment  transport 
rate  equation,  although  the  transport  rate  formulas  were  applied  in  a 
"2-line"  numerical  model  (Bailard  1985)  in  which  the  profile  is  characterized 
by  two  contours. 

59.  The  transport  model  of  Bailard  (1981)  generalized  Bagnold's 
steady-state  ene r ge t i c s- ba sed  stream  model  to  account  for  time-varying  flow 


and  an  arbitrary  sloping  bottom.  In  the  energetics  stream  approach,  it  is 
assumed  that  the  rate  of  sediment  transport  is  proportional  to  the  rate  of 
energy  dissipation  of  the  "stream, H  which  for  the  coastal  case  is  the  energy 
loss  caused  by  the  decay  of  broken  waves  and  by  bottom  friction. 

60.  Bailard's  formula  is  an  algebraic  expression  consisting  of  two 
terms  intended  to  describe  transport  contributions  from  bed  load  and 
suspended  load.  Empirical  efficiency  factors  enter  for  both  bed  and 
suspended  load;  these  factors  express  the  capability  of  the  stream  to  move 
the  sediment.  The  transport  equation  also  contains  the  sediment  fall 
velocity  (in  the  suspended  load  term),  the  local  beach  slope,  the  internal 
angle  of  friction  of  the  sediment  (for  bed  load),  and  higher  moments  of  the 
wave  orbital  velocity.  (In  the  general  case  in  which  longshore  transport  is 
included,  the  equation  contains  longshore  components  of  the  steady  and 
unsteady  wave  and  wave-induced  currents.)  The  wave  height  and  period  enter 
through  the  wave  orbital  velocities. 

61.  Bailard  (1982,  1983),  who  applied  his  transport  equation  to 
calculate  profile  change  observed  in  the  field,  reported  disappointing 
results,  finding  that  the  model  accounted  for  only  19  percent  of  the  observed 
variance  in  the  measured  beach  volume  change.  He  noted  that  actual  beach 
profile  change  is  difficult  to  calculate  directly  from  sediment  transport 
rates  because  the  net  transport  is  a  small  difference  resulting  from  two 
large  (onshore  and  offshore)  transport  rates.  The  required  moments  of  the 
wave  orbital  velocity  cannot  presently  be  calculated  from  existing  wave 
theory  with  any  confidence.  In  fact,  the  orbital  velocities  of  commonly  used 
small  amplitude  theory  are  symmetric  and  would  yield  no  net  transport  at  a 
given  location.  A  wave  theory  suitable  for  describing  surf  zone  waves  has 
not  yet  been  developed.  Bailard  used  empirical  moment  estimates  obtained 
from  current  meter  records.  These  instruments  and  records  may  not  be 
sufficiently  accurate  for  the  intended  purpose  (Aubrey  and  Trowbridge  1985). 

62.  In  summary,  the  transport  equation  of  Bailard  is  relatively 

soph i s t i ca t ed  and  ambitious  in  describing  many  of  the  physical  processes,  but 
its  applicability  has  not  been  demonstrated.  In  addition,  despite  its 
sophistication,  the  sediment  transport  equation  still  relies  on  empirical 
parameters  ("efficiencies,  moments,  and  friction  factor),  and  requires  know¬ 
ledge  of  quantities  which  are  difficult  to  calculate.  As  with  the  other 
models  introduced  above,  the  onshore  boundary  condition  has  not  been 


investigated  to  the  same  level  as  the  sediment  transport  rate  for  the 
subaqueous  profile.  We  also  note  that  the  energetics  approach  may  not  be 
applicable  seaward  of  the  wave  breaker  line. 

63.  At  present,  no  theoretical  model  of  beach  profile  change  or  dune 
erosion  exists  which  can  be  applied  for  engineering  purposes  now  or  in  the 
immediate  future.  Considerably  more  theoretical  and  numerical  groundwork 
must  be  laid.  Little  is  known  about  the  onshore  boundary  condition  required 
for  such  models.  In  addition,  sediment  transport  rate  predictive  equations 
would  be  best  evaluated  by  comparison  to  sediment  transport  rate  data  and  not 
solely  by  comparison  to  the  resultant  beach  profile  change. 


PART  IV:  EMPIRICAL  AND  PARAMETERIZED  MODELS 


64.  The  earliest  shoreline  retreat  model  was  proposed  by  Bruun  (1962) 
who  postulated  that  the  profile  would  rise  and  shift  landward,  without 
changing  shape,  in  response  to  a  rising  sea  level.  The  Netherlands  has 
always  been  concerned  with  storm  effects  and  dune  erosion,  and  Edelman  (1968, 
1972)  proposed  the  first  dune  erosion  model  based  on  a  simple  change  in 
profile  slope  and  a  balance  of  eroded  to  deposited  material.  Modified 
Edelman  approaches  have  been  used  by  Vallianos  (1975),  Tetra  Tech  (1983),  and 
by  Hill  and  Herchenroder  (1985). 

65.  Recent  research  efforts  have  resulted  in  two  models  which  are 
currently  being  used  for  engineering  studies.  Based  on  a  series  of  large- 
and  small-scale  laboratory  tank  tests,  limited  field  data,  and  theoretical 
work,  the  Dutch  have  improved  and  refined  the  Edelman  approach  (Vellinga 
1983b,  1986  and  van  de  Graaff  1983)  and  are  using  it  in  the  design  of 
artificial  dunes  in  The  Netherlands.  At  nearly  the  same  time,  a  numerical 
model  of  beach  and  dune  erosion  based  on  the  equilibrium  profile  equation  was 
developed  by  researchers  at  the  University  of  Delaware  and  the  University  of 
Florida  and  reported  in  a  series  of  theses  (Kriebel  1982  and  Moore  1982), 
papers  (Kriebel  and  Dean  1985a, b  and  Kriebel  1986),  and  reports  (Chiu  and 
Dean  1984  and  Kriebel  and  Dean  1985c).  The  model  has  been  widely  used  to 
estimate  dune  erosion  in  Florida  (Chiu  and  Dean  1984,  Kriebel  1984a, b,  and 
Kriebel  1986).  Both  of  these  models,  which  will  be  referred  to  as  the 
"Vellinga"  and  "Kriebel"  models,  have  been  shown  to  produce  reasonable  dune 
erosion  estimates.  A  third  model,  the  MSBWT  model  developed  by  Balsillie 
(1984c),  is  also  being  used  by  the  State  of  Florida  for  erosion  estimates. 

In  this  section,  these  three  models  will  be  described,  and  the  relative 
merits  of  each  will  be  discussed. 


The  Veiling 


Mode  1 


66.  Because  much  of  their  country  is  below  MSL,  the  Dutch  have  long 
considered  the  building  of  dikes,  dunes,  and  other  coastal  structures  of 
primary  importance.  Vellinga  (1983b,  1986)  presents  a  complete  discussion  of 
the  most  recent  Dutch  method.  A  probabilistic  approach  to  dune  design  using 
the  model  has  been  reported  by  van  de  Graaff  (1983,  1986).  Sargent  and 


moving  North  Sea  storm  surges  and  is  one  of  the  limiting  assumptions  of  the 
model  if  intended  for  general  use. 

68.  The  equation  is  applied  to  a  distance  200  to  250  m  seaward  of  the 
shoreline.  From  the  scale  relationships  used  in  the  model  studies,  Equation 
2  can  be  generalized  in  terms  of  Hos  and  fall  velocity  w  as 


where  x  ,  y  ,  and  Hos  are  in  meters  and  w  is  in  m/sec.  The  profile 
defined  by  Equation  3  terminates  offshore  at 


which,  substituting  into  Equation  3,  yields 

y  «  0.75  Hos  (5) 

69.  Constant  slopes  are  used  to  terminate  the  profile  at  the  shoreward 
end  (rri|  at  x  *  y  =  0)  and  seaward  end  (m2)  of  Equation  3  and  are  defined  as 

mi  =  -1 : 1  (6) 

m2  -  - 1 : 1 2 . 5  ( 7  ) 

Use  of  this  steep  poststorm  dune  face  mj  is  consistent  with  field 
observations,  and  m2  is  arbitrarily  fixed  in  agreement  with  model  tests. 

70.  Application  of  Equation  3  is  straightforward.  Based  on  the  wave 
height  and  fall  velocity,  a  profile  shape  is  computed.  The  origin  of  this 
line,  with  the  two  end  slopes  is  placed  at  the  surge  level  and  shitted 
horizontally  iandward  until  the  eroded  volume  equals  the  deposition.  In 


order  to  compensate  for  longer  duration  storms,  Vellinga  (1983b)  suggests 


that  the  erosion  volume  (above  the  surge  level)  be  increased  by  5  to  10 
percent  for  each  hour  (over  5)  that  the  storm  surge  is  within  1  m  of  the  peak 
surge  level  (with  a  maximum  increase  of  50  percent).  This  suggested  increase 
is  based  on  the  results  of  long  duration  model  tests  with  constant  surge 
leve 1 . 

71.  Note  that  as  a  result  of  the  simplicity  of  Equation  3,  the  post¬ 
storm  profile  shape  is  dependent  only  on  the  prestorm  profile,  wave  height, 
surge  level,  and  grain  size.  After  a  storm  passes,  only  a  larger  storm,  with 
a  higher  surge  level,  will  cause  additional  erosion.  The  model  is  based  on 
the  assumption  that  the  profile  adjusts  to  equilibrium  and  that  all  sediment 
movement  is  in  the  offshore  direction  (no  gradients  in  longshore  transport 
and  no  overwash).  The  model  does  not  account  for  the  formation  and  movement 
of  bars  and,  in  fact,  Vellinga  (1983b)  noted  that  the  model  overpredicted  the 
erosion  for  model  tests  with  offshore  bars. 

72.  Despite  these  limitations,  the  model  has  been  shown  to  produce 
reasonable  dune  erosion  estimates  (Vellinga  1983b,  1986,  and  Sargent  and 
Birkemeier  1985).  Based  on  all  available  model  and  field  data  (Figure  13), 
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Figure  13.  Predicted  volume  changes  above-surge  level  versus  measured 
laboratory  and  field  data  (Vellinga  1983b) 


Vellinga  (1983b)  specifies  the  standard  deviation  of  the  predicted,  above¬ 
surge  level,  eroded  volume  A„ 


as 


m  /m 


(8) 
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This  relationship  indicates  that  the  model  was  developed  and  is  most  useful 
for  large  storms  producing  above-surge  level  erosion  quantities  in  excess  of 
100  m^/m,  a  relatively  large  quantity  if  compared  to  published  data  on  eroded 
volumes  along  United  States'  coasts.  For  example,  using  above-MSL 
quantities,  only  one  of  the  storms  given  in  Table  1  caused  in  excess  of  100 
m^/m  of  erosion.  Similarly,  in  reports  by  Balsillie  (1985c,  1985d)  which 
document  the  effect  of  Hurricane  Elena  and  a  severe  northeaster  on  the  State 
of  Florida,  the  largest  single  profile  change,  above  MSL,  was  119  m^/m.  Most 
changes  were  considerably  less,  as  were  the  above-surge  changes. 

Properties  of  the  Vellinga  model 

73.  In  this  section  the  specific  properties  of  the  Vellinga  model,  as 
they  relate  to  the  requirements  of  a  dune  erosion  model  given  in  Part  III, 

are  addressed. 

a.  Deepwater  wave  height  is  included  both  in  the  equation  of  the 
predicted  profile  and  in  the  offshore  termination  point.  On 
the  basis  of  the  model  tests,  wave  steepness  had  limited 
effect  on  the  erosion  quantity.  The  model  is  applicable  for 
wave  steepnesses  between  0.02  and  0.04  (Vellinga  1983). 

b.  The  model  requires  specification  only  of  the  peak  surge  level. 
The  time-history  of  the  water  level  is  not  required.  Wave 
runup  and  setup  are  not  included.  The  duration  of  the  storm 
is  indirectly  simulated. 

c_.  The  shape  of  the  predicted  profile  is  dependent  on  the  fall 
velocity  of  the  median  grain  size;  the  native  slope  is  not 
directly  included. 

d^.  The  actual  prestorm  profile  shape  is  used.  This  profile 

allows  features  such  as  berms  and  bars  to  be  accounted  for  in 
the  balance  between  erosion  and  deposition.  However,  the 
required  specification  of  the  true  presto™  profile  is  a 
disadvantage  because  the  shape  of  the  profile  prior  to  the 
100-year  storm  will  not  be  known. 

e .  A  realistic  1:1  slope  is  assumed  i or  the  eroded  dune  from  the 
surge  level  up.  This  slope  is  consistent  with  field  measure¬ 
ments  of  nearly  vertical  beach  and  dune  scarping .  The  dune 
foot  is  specified  at  the  surge  level,  which  is  consistent  with 
data  from  Balsillie  (  1 '.-’3  5... )  .  Vellinga  (  1983b)  indicates  that 
field  data  have  shown  variation  of  the  dun-'  foot  from  0.5  to  - 


0.5  m  relative  to  the  surge  level.  The  model  does  not  predict 
berm  recovery. 


f . 


Because  the  model  is  based  on  an  assumed  equilibrium  of  the 
poststorm  shape,  the  rate  of  transport  is  not  directly 
specified.  Conservation  of  sediment  is  satisfied  by  requiring 
that  deposition  balance  erosion. 


8.- 


An  equilibrium  profile  is  computed  under  the  assumption  of  a 
5-hr  storm  surge  duration.  Longer  durations  are  treated 
indirectly.  This  approach  correctly  assumes  that  most  of  the 
upper  beach  and  dune  erosion  takes  place  during  the  time  of 
peak  surge. 


h. 


The  shape  of  the  poststorm  profile  has  been  verified  by  model 
tests  at  a  number  of  different  scales  and  for  field  profiles 
recorded  after  the  Dutch  storm  surge  of  1953  (Vellinga  1983b) 
Sargent  and  Birkemeier  (1985)  have  also  demonstrated 
applicability  of  the  model  for  modeling  erosion  caused  by 
minor  storms  in  the  United  States. 


74.  Although  the  Vellinga  model  lacks  a  detailed  description  of 
sediment  transport  and  is  limited  to  cross-shore  transport  on  sandy  beaches 
with  dunes,  as  will  be  shown  in  Part  V,  it  is  easy  to  apply  and  does  provide 
reasonable  estimates  of  dune  erosion.  Its  major  disadvantage  for  use  on 
United  States'  beaches  is  the  lack  of  a  duration  factor.  The  method  is 
expected  to  underestimate  dune  erosion  caused  by  long  duration  northeasters. 
In  addition,  a  poststorm  profile  shape  of  almost  fixed  form  limits  the 
applicability  of  the  procedure  since  it  prevents  the  method  from  working  on 
mildly  sloping  (flat)  beaches.  (It  was  found  to  be  impossible  to  balance  the 
erosion  and  deposition.) 


The  Kriebel  Model 


75.  The  Kriebel  model  is  based  on  the  equilibrium  profile  concept 
originated  by  Bruun  (1954)  and  further  investigated,  extended,  and  verified 
by  Dean  (1977),  Hughes  (1978),  Hughes  and  Chiu  (1981),  and  others.  The 
equilibrium  beach  profile  concept  is  discussed  in  Part  II.  The  fundamental 
assumption  of  the  model  is  that  the  profile  is  shaped  by  uniform  energy 
dissipation  per  unit  volume  of  breaking  and  broken  waves  in  the  surf  zone. 

By  this  principle,  the  model  is  limited  to  describing  the  evolution  of  the 
subaqueous  beach  profile  from  the  swash  zone  to  the  breaker  zone  (Figure  2). 
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Because  of  a  singularity  (described  below),  the  equilibrium  principle  cannot 
be  applied  at  the  shoreline;  the  equilibrium  concept-based  portion  of  the 
model  must  be  terminated  some  distance  seaward  of  the  water's  edge.  Ad  hoc 
geometrical-based  assumptions  are  employed  to  simulate  profile  change  in  the 


swash  zone  and  on  the  beach  face  and  dune.  Despite  these  apparent  severe 
limitations,  the  Kriebel  model  has  been  demonstrated  to  yield  qualitatively 
and  quantitatively  reasonable  results  (Kriebel  1982,  1986). 

Governing  equations 

76.  The  following  discussion  is  based  on  the  work  of  Dean  (1977), 
Kriebel  (1982),  Moore  (1982),  and  Kriebel  and  Dean  (  1985a, b).  The  two 
reports  by  Kriebel  (1984a, b)  are  recommended  for  a  more  complete  discussion. 
The  fundamental  assumption  is  that  the  beach  profile  is  well  approximated  by 
the  following  relationship: 

h(x)  =  Ax*3  (9) 

which  is  the  general  power  law  form  of  Equation  1.  As  discussed  in  Part  II, 
the  most  appropriate  value  of  the  scale  coefficient  b  is  2/3,  which  can  be 
explained  by  uniform  energy  dissipation  (Dean  1977).  However,  use  of  a 
variable  exponent,  although  highly  empirical,  results  in  better  fits  to  some 
profiles  (Balsillie  1982).  The  best  fit  scale  coefficients  determined  in 
this  manner  do  not  deviate  greatly  from  the  value  of  2/3.  The  shape 
parameter  A  was  empirically  related  to  the  mean  grain  size  (Figure  8)  by 
Moore  (1982).  Moore  recommended  that  A  be  represented  by  a  smooth  line 
drawn  through  somewhat  scattered  results  to  give  a  monotonica 1 ly  increasing 
function  of  grain  size.  This  relationship  expresses  the  well-known  result 
that  coarse-grained  beaches  tend  to  be  steeper  than  fine-grained  beaches. 

77.  The  basic  equation  for  the  cross-shore  transpo-'^  rate  Qc  per 
unit  beach  width  is  assumed  to  be  given  by 


where 


QC  “  k(D  Dan ) 


=  empirical  coefficient 

*  energy  dissipation  per  unit  volume  (assuming  spilling  waves) 

=  energy  dissipation  corresponding  to  the  equilibrium 
profile  for  a  beach  of  given  grain  size 


Operationally,  the  spilling  wave  assumption  means  that  the  wave  height  of 
breaking  and  broken  waves  is  proportional  to  the  local  water  depth,  using  a 
fixed  constant  of  proportionality.  With  the  x-axis  directed  offshore  as 
shown  in  Figure  7,  offshore  transport  has  a  positive  sign,  and  onshore 
transport  has  a  negative  sign.  After  manipulations,  the  dissipation  D  is 
ultimately  expressed  in  the  form 

D  -  const  ^h1/2  (11) 


in  which  const  is  the  product  of  known  constant  factors.  The  empirical 
coefficient  k  was  found  to  have  the  value  2.2  x  10“^  m^/N  (Moore  1982) 
based  on  the  results  of  profile  change  found  in  the  large-scale  flume 
experiments  of  Saville  (1957).*  From  Equations  10  and  11,  it  is  seen  that 
if  the  depth  is  greater  than  the  equilibrium  depth  at  a  given  location,  sand 
will  move  offshore,  i.e.,  erosion  is  associated  with  the  higher  water  levels 
which  would  occur  during  a  storm. 

78.  It  is  noted  that  wave  height  and  wave  period  do  not  explicitly 
enter  in  the  transport  rate  equation.  The  water  level  (the  depth  h)  is  seen 
to  be  the  main  external  force  determining  the  cross-shore  transport  rate. 

The  wave  height  is  indirectly  included  in  the  transport  equation  since  it  is 
used  to  determine  the  location  of  the  breaker  depth  or  width  of  the  surf  zone 
over  which  the  model  acts. 

79.  Beach  profile  change  is  put  into  time-dependent  form  by  inserting 
the  transport  rate  predictive  formula,  Equation  10,  into  the  beach  material 
(sand)  continuity  equation  expressed  in  the  form 


_  _  ?2c 

dt  dh 


(12) 


in  which  t  is  time.  Equation  12  can  be  numerically  solved  if  the  surge 
hydrograph,  initial  profile,  and  offshore  wave  height  and  period  are  given 
The  hydrograph  describes  the  water  level  as  a  function  of  time  and  must  be 
either  predicted  by  another  numerical  model  or  be  specified  from  measure¬ 
ments  . 

*  Kriebel  (1986)  recalibrated  his  revised  model  using  data  from 
Saville  (1957)  and  determined  that  k  should  equal  8.7  *  10”^  m^/N. 
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80.  If  Equation  10  is  inserted  into  Equation  12,  one  term  results  with 
an  exponent  of  the  depth  h  as  -1/2.  This  quantity  is  singular  at  the 
shoreline  where  h  «  0.  Therefore,  in  practice,  the  equilibrium  profile- 
based  portion  of  the  model  is  terminated  at  some  finite  water  depth  and 
connected  to  another  model  which  operates  landward  of  the  connection  point. 
The  Kriebel  model  uses  simple  geometrical  and  volume  conserving  procedures  in 
this  landward  section.  In  particular,  a  linear  beach  face  slope  is  smoothly 
connected  with  the  equilibrium  profile  at  the  depth  where  the  beach  face 
slope  equals  the  slope  of  the  equilibrium  profile  for  the  given  water  level 
conditions.  In  the  version  of  the  model  evaluated  in  Part  V  (Kriebel 

1984a, b),  two  options  are  permitted  for  simulating  dune  erosion:  one  allows 
representation  of  a  wide,  flat  berm  and  the  other  a  sloping  beach  that 
directly  connects  to  a  dune.  Kriebel  (1986)  includes  a  more  realistic  wave 
runup  limit,  the  formation  of  an  erosional  scarp,  and  a  more  realistic, 
nearly  vertical  eroded  dune  face. 

Properties  of  the  Kriebel  model 

81.  The  Kriebel  model  satisfies  a  number  of  the  criteria  developed  in 
Part  III  for  judging  the  quality  of  a  dune  erosion  model.  A  critical 
discussion  will  now  be  given  based  on  general  properties  of  the  model. 

a.  Wave  height  and  period  do  not  explicitly  appear  in  the  model. 
The  wave  height  is  used  only  in  the  determination  of  the  wave 
breaking  point,  hence  the  width  of  the  surf  zone.  These 
important  parameters  effectively  disappear  because  of  the 
basic  assumptions  of  the  model,  i.e.,  equilibrium  profile, 
shallow-water  approximations  to  small  amplitude  wave  theory, 
and  the  spilling  wave  breaker  assumption.  Since  wave 
steepness  or  another  criterion  is  not  used  to  specify  whether 
erosion  or  accretion  will  occur,  the  model  may  predict  erosion 
when  actually  accretion  or  recovery  would  occur  for  the  given 
wave  and  beach  conditions. 

b.  The  model  requires  specification  of  the  water  level  through 
time  which  it  uses  in  a  time-stepping  manner.  Both  storm 
duration  and  water  level  enter  through  the  characteristics  of 
the  hydrograph. 

c.  The  model  accounts  for  grain  size  or  fall  velocity  through  the 
scale  parameter  A.  Values  of  A  were  determined  empirically 
and  thus  can  be  used  with  relative  confidence. 

d.  The  initial  profile  is  not  required,  although  it  can  be  used 
if  available.  This  is  both  a  potential  advantage  and 
disadvantage.  On  the  one  hand  ,  an  initial  profile  will 


probably  not  be  accurately  known,  so  that  use  of  a  schematic 
profile  is  convenient.  On  the  other  hand,  beach  erosion  is 
expected  to  be  sensitive  to  the  initial  profile  configuration, 
including  the  presence  of  a  bar. 

e The  onshore  boundary  condition  is  based  on  simple  geometrical 
considerations.  It  does  not  accurately  reproduce  poststorm 
profile  shapes.  A  more  sophisticated  onshore  boundary 
condition  is  desirable. 

f_.  A  highly  schematic  but  very  reasonable  expression  is  used  to 
model  the  cross-shore  sediment  transport  rate.  Although  not 
based  on  fundamental  physical  principles,  it  appears  to  be 
well  suited  for  engineering  applications.  The  sediment 
continuity  equation  is  used  to  calculate  profile  change  and 
dune  erosion.  It  is  limited  to  cross-shore  transport,  and  the 
only  major  features  modeled  are  dune  and  berm  erosion  and 
change  of  an  assumed  equilibrium-shaped  profile. 

£.  The  model  time-steps  through  the  storm  hydrograph  toward  an 
equilibrium  form.  This  time-dependent  capability  is  an 
important  property  of  the  model  since  it  allows  both  long  and 
short  duration  storms  to  be  modeled. 

h.  The  model  has  been  demonstrated  to  give  reasonable  results  in 
a  limited  number  of  comparisons  to  measured  dune  erosion 
volumes.  A  more  detailed  discussion  of  this  point  is  given  in 
Part  V. 


31  .  The  Kriebel  model  is  judged  to  provide  a  reasonable  means  to 
estimate  dune  erosion  which  is  compatible  with  the  expected  quality  of  the 
input  data.  The  model  requires  as  input  the  beach  grain  size;  an  initial 
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profile  if  available  or,  instead,  an  assumed  schematic  profile;  the  surge 
hyd  r -■g -a;  h  ;  and  the  deepwater  wave  height.  The  model  does  not  allow  for  bar 
format  ior. ;  therefore,  the  recovery  process  is  not  well  described.  Similarly, 
:;ed  i  mer;'  transport  seaward  of  the  breaker  line  is  not  taken  into  account. 

3'.  Although  not  discussed  in  detail  here,  the  Kriebel  model  could  be 
modified  t o  allow  for  inundation.  In  principle,  the  model  could  also  be 
ext  ended  *o  describe  erosion  of  seawall-backed  dunes.  Such  a  model  is 
-ur :  e:.t  !  y  undergoing  testing  at  CERC  and  is  found  to  perform  favorably. 


r.+ 


34.  The  MSBWT  model  (Balsillie  1984c)  resulted  from  a  series  of 
investigations  of  breaking  wave  properties  (Balsillie  1984a, b)  in  order  to 
estimate  the  height  of  storm  waves  crossing  low-lying  inland  areas. 

Recently,  the  erosion  process  has  been  added  to  the  model,  and  it  is  now  used 
to  est  in. ate  storm  effects  on  nonflooded,  flooded,  and  breached  profile  types 
(Balsillie  1935b). 

G  o  v  e  mi  n  g  e  qua  t  ions 

85.  The  MSBWT  model  is  based  on  a  combination  of  the  physics  of  surf 


zone  wave  dynamics  and  statistical  estimates  of  unknowns  such  as  erosion 
quant  it  i e  s  and  bar /tn  ugh  shapes.  The  basic  assumption  of  the  wave 
t : «r:sf or mat  ion  model  is  that  energy  in  the  surf  zone  is  not  uniformly 
dissipated  but  is  dissipated  through  the  breaking  and  reformation  of  incident 
waves  (a  ir  're  realistic  assumption).  Breaking  waves  are  required  to  be 
plunging  breakers.  This  is  more  appropriate  than  the  spilling  wave 
assumption,  since  waves  typically  break  by  plunging  on  most  coasts  of  the 
United  States. 

86.  The  model  in  based  on  the  assumption  of  an  initial  equilibrium 
profile  shape  for  the  offshore  and  uses  the  actual  beach  shape.  Beginning 
with  the  surge  height  (including  wave  setup),  a  design  wave  period,  and  a 
20-fr  (6-m)  wave,  waves  in  1-ft  (0.3  m)  high  increments  are  propagated  shore¬ 
ward.  At  computed  plunge  points,  an  offshore  bar  and  trough  are  formed,  the 
shapes  of  which  are  determined  by  the  statistical  analysis  of  bar/trough  data 
for  the  depth  of  the  particular  breaker. 

37 .  Since  the  bars  form  from  the  convergence  of  sediment  both  from 
offshore  and  onshore,  a  net  deficit  of  material  exists  at  the  trough  of  the 
innermost  tar.  The  depth  of  this  trough,  plus  an  additional  depth  based  on 
the  depth  of  bed  liquefaction,  is  used  to  generate  a  poststorm  profile  shape 
extending  up  to  the  surge  level.  Similar  to  the  Vellinga  model,  a  1:1  slope 
is  assumed  for  the  face  of  the  poststorm  dune  (above  the  surge  level). 

88.  The  eroded  volume  is  separately  computed  based  on  the  following 
dimensional  (metric)  equation  (Balsillie,  1985a,  1986): 
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where 

Qe  aVg  =  average  erosion  for  a  particular  area,  m^/m 

g  =  acceleration  of  gravity,  m/sec^ 

tr  =  rise  time  (the  time  required  to  reach  the  peak  surge 
level ) ,  sec 

S  =  peak  storm  surge,  m 

The  data  from  the  10  cases  used  to  determine  this  equation  show  a  remarkable 
linear  fit  (Figure  14).  This  equation  is  interesting  because  it  relates  the 
final  erosion  quantity  only  to  the  surge  level  and  the  rise  time  rather 


Figure  14.  Relationship  between  measured  average  net  erosion 
Qe  aVg  and  the  factor  (g^^  tr  S^)  (Balsillie  1986) 

than  to  grain  size,  beach  slope,  offshore  features,  or  wave  height.  Unlike 
the  Kriebel  and  Vellinga  models,  since  the  profile  does  not  adjust  to 
equilibrium  based  on  a  balance  between  erosion  and  deposition,  successive 
similar  storms  will  produce  the  same  amount  of  erosion.  As  in  the  Vellinga 
model,  the  peak  surge  is  most  important.  The  duration  is  incorporated 
through  use  of  the  rise  time,  which  correctly  reduces  the  erosion  potential 
of  fast  moving  storms. 

89.  Since  for  a  particular  purpose  the  model  of  Balsillie  (1985a)  was 


v-'-.v-v-v 


v*  *  *  -  j.  -  x  ■  •  v*  *•  hV.  '  A  A 


'll  IE 


d  to  predict  the  92.5  percentile  of  erosion,  the  following  equat i 
sed  : 


Qe  max 


nee  the  erosion  quantity  is  determined  from  Equation  13  (or  Equation  19), 
he  poststorm  shape  is  moved  along  the  surge  level  until  the  eroded  area 
qua  Is  the  estimate. 

90.  E  e  cause  the  MSBWT  model  is  designed  to  predict  the  average  (or 

aximurni  erosion  for  a  particular  area,  it  was  verified  (Balsillie  1985b)  bv 
onqurirg  actual  and  predicted  poststorm  profile  shapes  to  32  cases  where 
“as. re i  changes  were  within  i.5  rn^/m  of  the  actual  Qe  av  In  addition, 

a mage  to  three  piers  i  two  ,n  Florida,  one  in  California)  was  correlated  to 
he  maxim  ur,  wave  crest  elevation  predicted  by  the  model, 
r  c'er-tie.-,  of  the  HSFWT 

91.  This  section  relates  the  properties  of  the  MSBWT  model  to  the  list 
f  requirements  for  dune  erosion  models  stated  in  Part  III. 

a.  The  deepwater  wave  height  and  period  entered  are  required  by 
the  wave  transformation  model  in  order  to  determine  tne  depth 
of  the  innermost  trough.  Wave  steepness  is  included  through 
use  of  a  surf  parameter  defined  as  the  wave  steepness  divided 
by  the  square  root  of  the  beach  slope  (Balsillie  1984c). 

b.  The  model  requires  input  of  the  maximum  surge  level  above  the 
still-water  line,  including  the  wave  setup. 

r_.  The  effect  of  grain  size  is  included  only  as  it  affects  the 
offshore  slope  and  the  predicted  equilibrium  post  storm  shape. 

d_.  The  actual  beach  profile  shape  is  fitted  to  an  equilibrium 
offshore  shape.  The  true  shape  of  the  offshore  is  not 
required  . 

e.  The  onshore  boundary  is  dependent  on  the  type  of  profile 
(nonflooded,  flooded,  or  breached).  Each  is  treated 
differently.  On  nonflooded  profiles,  a  1:1  slope  is  used  for 
the  eroded  dune  above  the  surge  level. 

f.  Cross-shore  transport  is  modeled  through  the  volumetric 
requirements  of  bar  formation.  The  model  has  provisions  for 
accounting  both  dune  erosion  and  over wash .  Volume  losses  ate 


determined  by  an  independent  equation  so  that  sand  volume  is 
not  conserved. 

g_.  The  wave  transformation  model  is  not  time  dependent;  it 
assumes  that  the  offshore  bar-trough  bathymetry  is  in 
equilibrium.  The  erosion  model  includes  time  dependency 
through  use  of  the  rise  time  of  the  surge  level. 

h.  The  model  has  been  tested  for  a  range  of  Florida  and  gulf 

coast  profile  changes  and  found  to  produce  reasonable  results. 
The  volumetric  relationship  (Equations  13  and  14)  is  based  on 
a  limited  number  of  data  points  but  from  a  wide  range  of 
Florida  and  non-Florida  beaches. 

92.  The  MSBWT  model  has  evolved  from  a  surf  zone  wave  transformation 
mode  1  to  a  combined  wave  and  beach  profile  erosion  model.  Because  the  model 
includes  a  description  of  many  of  the  important  processes  associated  with 
dune  erosion,  further  testing  and  evaluation  should  be  made,  particularly 
with  respect  to  non-Florida,  non-gulf  coast  beaches.  Several  considerations 
preclude  use  of  the  MSBWT  at  present.  Although  numerous  reports  document 
various  aspects  of  the  model,  a  report  by  Balsillie  (in  preparation) 
detailing  the  present  model  in  use  by  the  State  of  Florida  has  not  been 
published.  Though  the  code  for  the  model  has  been  published  (Balsillie 
1984-'  i,  it  is  written  in  the  Applied  Programming  Language,  a  computational 
programming  language  which  is  seldom  encounteied  in  coastal  research. 
Conversion  of  the  complex  code  to  FORTRAN  or  PASCAL  would  be  time  consuming 
ar.d  costly.  Finally,  the  model  and  the  equations  and  assumptions  on  which  it 
is  based  must  pass  the  scrutiny  of  independent  testing  and  ver i f icat ion . 


PART  V:  EVALUATION  OF  THE  VELLTNGA  AND  KRIEBEL  MODELS 


93.  Both  the  Vellinga  and  Kriebel  models  are  highly  empirical  in  their 
conception  and  formulation.  However,  empirical  modeling  techniques  have  been 
proven  to  be  capable  of  predicting  certain  extremely  complex  natural  phenom¬ 
ena,  such  as  dune  erosion.  These  predictions  can  be  reliable  if 

a.  The  modeling  parameters  have  been  carefully  selected  for  a 
particular  site . 

b.  The  model  results  have  been  properly  interpreted  in  a  manner 
consistent  with  the  assumptions  and  simplifications  inherent 
in  the  model  formulation. 

c_.  The  modeled  area  and  storm  event  do  not  significantly  deviate 
from  the  conditions  for  wh^ch  the  model  was  formulated  and 
tested. 

94.  Since  a  certain  amount  of  interpretative  skill  and  engineering 
judgment  are  required  for  use  of  both  the  Vellinga  model  and  the  Kriebel 
model,  a  quantitative  comparison  of  the  two  models  was  made  to  demonstrate 
their  capabilities.  Sargent  and  Birkemeier  (1935)  used  a  number  of  indi¬ 
vidual  beach  and  offshore  profiles  for  each  of  four  storms  to  demonstrate  the 
accuracy  of  the  Vellinga  model.  In  this  section,  the  Kriebel  model  is 
applied  to  these  same  data  and  compared  both  to  the  Vellinga  model  results 
and  to  measured  changes.  The  comparison  is  based  on  only  14  points.  A  more 
rigorous  evaluation  with  a  wider  variety  of  data,  though  desirable,  was 
beyond  the  scope  of  this  feasibility  study. 

Field  Data 

95.  The  storm  data  used  by  Sargent  and  Birkemeier  (1985)  are  summar¬ 
ized  in  Table  2.  Deepwater  wave  data  were  either  obtained  from  3-hr  wave 
hindcasts  (Jensen  1983)  in  9  m  of  water  ( Westhampton ) ,  from  gage  measurements 
in  comparable  water  depths  (Duck),  or  estimated  from  visual  breaking  wave 
observations  (Long  Beach  Island).  Surge  heights  (without  wave  setup  or 
runup)  were  computed  from  water  level  measurements  collected  by  the  nearest 
available  sensors.  Sediment  fall  velocities  were  determined  from  repre¬ 
sentative  sediment  grain  size  data  from  each  beach. 

96.  The  severity  of  each  of  the  four  storm  events  used  in  the  compar- 
i  ■  >r.  was  quantified  by  both  pre-  and  poststorm  profile  surveys,  including 


Table  2 


Storm  Events 


locatlon/Date 

buratlon* 

hr 

Surge  Level 
m  above  MSL 

Wave  Ht 

HI 

Fall  Vel 
m/sec 

Number  of 
Profiles 

Westhampton,  N.V. 

3  Feb  1972 

6.0 

1.3 

4.3 

0.0450 

3 

Westhampton,  N.Y. 

19  Feb  1972 

7.0 

1.5 

5.5 

0.0450 

3 

Long  Beach  Is,  N.J. 
19  Dec  1977 

9.0 

1.4 

3.0 

0.0450 

5 

Duck,  N.C. 

14  Nov  1981 

15.5 

1.6 

4.0 

0.0500 

3 

*  This  Is  duration  as  defined  by  the  Velllnga  model,  the  elapsed  time  the 
surge  level  was  within  1  m  of  the  peak  level. 


either  actual  (Duck)  or  approximate  offshore  data,  at  several  locations. 

Note  that  these  data  were  originally  selected  to  satisfy  the  requirements  of 
the  Vellinga  model  with  respect  to  dune  height,  width,  and  slope.  This  data 
set  was  further  analyzed  in  order  to  derive  the  input  conditions  required  for 
use  of  the  Kriebel  model.  This  analysis  involved  approximating  the  modeled 
area  by  defining  a  constant  dune  and  berm  height,  dune  and  berm  face  slopes, 
and  berm  width.  The  offshore  profile  was  approximated  as  monotonically 
increasing  in  depth  according  to  the  equilibrium  profile  concept. 

97.  The  resulting  approximations  for  the  geometry  of  each  of  the  sur¬ 
veyed  profiles  which  were  used  as  input  for  the  Kriebel  model  are  presented 
in  Table  3.  The  value  of  the  shape  coefficient  A  for  each  of  the  profiles 
was  selected  to  best  represent  the  prestorm  offshore  profile.  Since  the 
majority  of  the  profiles  are  characterized  by  bar  formations,  the  resulting 
equilibrium  profile  represents  a  smoothing  of  the  existing  offshore 
bathymetry . 

98.  An  additional  parameter  required  by  the  Kriebel  model  is  the  surge 
hydrograph,  the  duration  of  which  differs  from  that  defined  by  the  Vellinga 
model.  The  period  used  for  this  application  represents  the  duration  of  the 
entire  storm  surge  event;  whereas  the  duration  used  for  the  Vellinga  model 
only  represents  the  maximum  peak.  Therefore,  the  durations  shown  in  Table  3 
are  much  longer  than  those  of  Table  2. 
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Table  3 

Geometric  Input  Data 
Kriebel  Model 


Dune  Berm 


Locatlon/Oate 

Profile  Height 
Number  m 

Slope 

Width 

m 

Height 

m 

Slope 

ml'3 

Duration 

hr 

Westhampton 

3 

5.8 

0.174 

4.0 

3.4 

0.083 

0.124 

60 

3  Feb  1972 

4 

4.7 

0.070 

3.7 

2.0 

0.065 

0.114 

60 

5 

4.7 

0.098 

18.0 

2.0 

0.105 

0.075 

60 

Westhampton 

3 

5.6 

0.138 

0.0 

2.4 

0.066 

0.108 

60 

19  Feb  1972 

4 

4.6 

0.157 

0.0 

1.8 

0.040 

0.117 

60 

5 

8.4 

0.151 

10.0 

2.7 

0.051 

0.128 

60 

Long  Beach 

14 

5.5 

0.340 

0.0 

2.6 

0.086 

0.114 

96 

Island 

15 

4.3 

0.175 

0.0 

2.3 

0.078 

0.113 

96 

19  Dec  1977 

16 

5.3 

0.156 

4.9 

2.0 

0.087 

0.114 

96 

17 

5.0 

0.333 

0.0 

3.1 

0.081 

0.118 

96 

18 

5.3 

0.130 

2.4 

2.0 

0.109 

0.124 

96 

Duck 

186 

7.5 

0.210 

0.0 

1.7 

0.046 

0.105 

72 

14  Nov  1981 

188 

6.6 

0.203 

0.0 

2.4 

0.079 

0.101 

72 

190 

6.4 

0.500 

0.0 

4.0 

0.107 

0.101 

72 

Model  Comparison 

99.  Above-MSL  and  above-surge  level  erosion  quantities  predicted  by 
the  Vellinga  and  Kriebel  models  for  each  profile  line  are  given  in  Table  4. 
Note  that  as  recommended  by  Vellinga,  his  prediction  has  been  adjusted  for 
the  duration  of  the  storm,  which  in  most  cases  improves  the  result.  Although 
there  is  considerable  scatter  in  the  data  (Figure  15),  the  volumetric  esti¬ 
mates  are  deemed  reasonable.  The  worst  case  for  both  models,  Profile  18  on 
Long  Beach  Island,  is  an  isolated  case  which  probably  resulted  from  the 
location  of  this  profile  near  a  large  terminal  groin.  For  the  14  cases,  the 
Vellinga  model  tends  to  equally  over-  and  underpredict  the  above-MSL  changes 
but  underpredicts  all  but  two  of  the  above-surge  level  quantities.  The 
Kriebel  model  tends  to  overpredict  both  quantities. 

100.  Poststorm  profile  shapes  are  not  as  accurately  reproduced. 
Cross-section  plots  for  all  the  cases  are  shown  in  Figure  16.  Note  that 
because  only  the  poststorm  profile  shapes  from  the  Kriebel  model  are  shown  in 
Figure  16,  they  appear  misleading  with  respect  to  the  conservation  of  sand. 
Figures  17a  and  b  illustrate  the  schematic  initial  and  predicted  profile 
shapes  used  with  the  Kriebel  model  for  both  the  best  and  worst  fits  given  in 


Table  4. 

101.  The  data  presented  in  Table  4  indicate  that  both  models  are 


Table  4 


Measured  and  Predicted  Erosion  Quantities 


Erosion  Above  F6L 


Vellinga  Model 


Kriebel  Model 


Erosion  Above  Peak  Serge  Level 
Veil  Inga  Model  Kriebel  Model 


Location/ 

Date 

Profile 

Nurber 

Actual 

Predicted  Adjusted” 

Percent 

Actual 

Percent 
Predicted  Actual 

Actual 

Predicted  Adjusted” 

Percent 

Actual 

Percent 
Predicted  Actual 

westhanpton 

3 

52.8 

38.5 

41.1 

78 

29.5 

56 

18.2 

14.2 

15.6 

86 

20.3 

112 

3  Feb  1972 

4 

69.4 

11.6 

12.2 

18 

23.0 

33 

42.6 

2.3 

2.5 

6 

10.2 

24 

5 

51.3 

37.9 

39.8 

78 

32.3 

63 

33.8 

15.6 

16.4 

49 

19.2 

57 

westhanpton 

3 

26.5 

43  0 

47.3 

178 

33.9 

128 

15.3 

13.7 

15.1 

99 

24.8 

162 

19  Feb  1972 

4 

35.2 

10.5 

11.6 

33 

19.1 

54 

8.5 

0.2 

0.2 

2 

9.2 

106 

5 

37.1 

20.5 

22.6 

61 

21.1 

67 

17.1 

2.6 

2.9 

17 

7.7 

45 

Long  Beach 

14 

30.6 

26.2 

31.4 

103 

38.2 

125 

17.6 

9.1 

10.9 

62 

22.6 

128 

Island 

15 

26.2 

27.6 

33.1 

126 

42.5 

162 

18.2 

8.9 

10.7 

59 

30.8 

169 

19  Dec  1977 

16 

13.8 

25.8 

31.0 

224 

35.2 

255 

8.0 

5.9 

7.1 

89 

16.3 

204 

17 

34.5 

27.4 

32.9 

95 

43.4 

126 

20.0 

7.7 

9.2 

46 

33.6 

168 

18 

6.0 

20.8 

25.0 

416 

41.7 

695 

0.0 

7.3 

8.5 

- 

29.0 

- 

Duck 

186 

31.7 

20.6 

30.9 

97 

32.7 

103 

10.1 

3.0 

4.5 

45 

23.1 

229 

14  Nov  1961 

188 

41.3 

34.5 

51.8 

125 

48.1 

112 

18.5 

13.6 

20.4 

110 

37.8 

204 

190 

37.4 

44.0 

66.0 

176 

53.7 

143 

19.3 

20.3 

30.5 

158 

44.9 

233 
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event,  since  entirely  different  assumptions  were  made  in  their  development. 
An  obvious  example  is  the  difference  in  surge  duration;  nevertheless,  one 
type  of  comparison  of  the  models  can  still  be  made.  This  comparison  is 
accomplished  by  examining  the  sensitivity  of  both  models  using  the  percent 
changes  in  predicted  volume  of  erosion  which  occur  as  a  result  of  similar 
percentage  changes  in  equivalent  input  parameters.  Table  5  presents  several 
of  these  comparisons  which  were  computed  from  the  data  presented  in  Figures 
18  through  22  demonstrating  that  both  models  are  stable  for  the  parameter 
range  tested.  The  Vellinga  model  is  most  sensitive  to  the  surge  level  used 
and  the  grain  size.  The  Kriebel  model  is  less  sensitive  to  changes  in  most 
parameters . 


Table  5 

Sensitivity  Testing 


Figure 

Number 

Constants* 

Varlabl 

Change  In  Variable 

Percent  Increase  In  Erosion 
Above  Surge  Level 

a  From 

To 

Percent 

Vellinga  Kriebel 

18 

S 

.  2.5 

X 

O 

6.0 

°50 

0.25 

0.40 

60 

130 

18 

°so 

0.35 

5 

2.5 

Hos 

4.0 

8.0 

100 

40 

19 

Hos 

6.0 

U50 

0.35 

s 

2.0 

4.0 

100 

190 

19 

S 

3.0 

°50 

0.35 

Hos 

4.0 

8.0 

100 

80 

20 

Dur 

60 

S 

3.0 

A 

0.21 

o 

o 

0.73 

Hos 

4.0 

8.0 

100 

10 

21 

Hos 

6.0 

Dur 

60 

S 

3.0 

A 

0.079 

0.214 

171 

-70 

21 

Hos 

6.0 

Dur 

60 

A 

0.15 

o 

U"» 

o 

0.33 

s 

2.0 

4.0 

100 

30 

21 

Hos 

6.0 

Our 

60 

A 

0.21 

o 

O 

0.73 

s 

2.0 

4.0 

100 

110 

22 

Hos 

6.0 

5 

3.0 

A 

0.15 

°so 

0.33 

Dur 

60 

108 

80 

50 

*  Units:  HQS  and  S  In  qi,  duration  (Dur)  In  hr,  050 
A  In 


In  wn,  and 


Other  Evaluations 


103.  Recently,  both  Kriebel  (1986)  and  Vellinga  (1986)  presented 
additional  evaluation  data  for  their  models.  Kriebel,  using  an  improved  and 
recalibrated  version  of  the  model  tested  here,  evaluated  20  poststorm 
Hurricane  Eloise  profiles  from  Walton  County,  Florida,  which  were  deemed 
representative  of  the  80  to  90  eroded  profile  lines  available.  Maximum 
measured  erosion  above  MSL  was  approximately  60  m^/m.  All  20  cases  were 
within  a  40-percent  margin  of  error,  with  17  cases  lying  within  25  percent. 
Dune  recession  tended  to  be  underpredicted. 

104.  In  his  dissertation,  Vellinga  (1986)  presents  comparisons  for 
three  different  events,  including  one  profile  from  Hurricane  Eloise.  The 
largest  storm,  the  1953  storm  surge,  caused  an  average  of  90  m^/m  of  erosion 
and  formed  the  basis  for  much  of  the  development  of  the  Vellinga  and  earlier 
Dutch  models.  The  Vellinga  model  estimated  an  average  loss  of  97  m^/m  for 
this  event.  The  comparison  containing  the  most  data  (58  cases)  documented 
the  1976  storm  surge  which  produced  an  average  loss  of  32  m^/m.  This  overall 
loss  is  well  estimated  by  the  model,  though  individual  profile  changes  show 
considerable  scatter  with  differences  up  to  400  percent.  Good  agreement  was 
also  obtained  for  the  one  profile  from  Hurricane  Eloise. 

105.  Chiu  and  Dean  (1986)  evaluated  two  Florida  profiles  eroded  by 
hurricanes  using  a  highly  simplified  version  of  the  Kriebel  model  (Chiu  and 
Dean  1984)  and  found  that,  even  with  a  variability  factor  of  2.5,  the  pre¬ 
dicted  horizontal  recession  of  the  dune  face  was  underestimated. 

Cone lusions 

106.  Based  on  the  comparison  in  the  above  section,  both  the  Vellinga 
and  Kriebel  models  have  been  shewn  to  produce  reasonable  estimates  of  dune 
erosion.  The  models  are  stable  with  respect  to  the  input  parameters  if  they 
lie  within  reasonable  physical  limits.  The  Vellinga  model  requires  a  detail 
ed  knowledge  of  the  initial  beach  and  offshore  profile  and  the  storm  surge 
level  in  order  to  predict  storm  erosion  volumes  and  the  associated  poststorm 
profiles.  It  can  be  seen  from  the  governing  equation  of  the  Vellinga  model 
that  both  the  predicted  profile  shape  and  offshore  termination  point  are 
highly  dependent  on  the  estimated  deepwater  significant  wave  height. 


Moreover,  the  method  is  expected  to  work  best  with  profiles  similar  in  slope 
to  those  on  the  Dutch  coast.  As  pointed  out  by  Sargent  and  Birkemeier 
(1985),  for  much  flatter  slopes  the  predicted  profile  may  be  too  steep  to 
balance  the  erosion  and  deposition.  Finally,  an  accurate  representation  of 
the  offshore  profile  shape  through  the  storm  surf  zone  is  required.  The  more 
different  the  prestorm  and  predicted  poststorm  shapes  are,  the  higher  the 
erosion.  Offshore  survey  data,  particularly  through  the  surf  zone,  are 
difficult  and  costly  to  obtain,  and  the  dynamic  nature  of  this  zone  dimin¬ 
ishes  the  value  of  a  single  survey  for  long-term  estimates  by  FEMA.  Use  of 
the  Vellinga  model  is  not  recommended  except  possibly  to  supplement  the 
Kriebel  model  for  conditions  suited  to  its  range  of  applicability.  It  could 
be  particularly  useful  in  situations  where  the  prestorm  offshore  profile 
shape  is  known  and  the  surge  duration  is  not. 

107.  The  Kriebel  model  is  a  more  sophisticated  and  generalized  dune 
erosion  predictive  technique  than  the  empirically  based  Vellinga  model, 
because  it  is  based  on  a  predictive  transport  rate  formula  and  has  a  time 
dependency.  The  Kriebel  model  requires  only  a  schematized  prestorm  profile 
(which  can  be  obtained  from  measurements  or  calculated  within  the  framework 
of  the  model),  storm  surge  hydrograph,  and  the  deepwater  wave  height.  The 
capability  to  use  a  generalized  profile  shape  is  an  advantage  compatible  with 
FEMA's  standardized  calculation  procedures  because  the  model  is  relatively 
insensitive  to  small  changes  in  the  beach  profile  that  continuously  occur 
through  time.  In  conclusion,  based  primarily  on  this  ability  to  use  a 
schematic  prestorm  profile  and  because  it  is  relatively  insensitive  to  minor 
variations  in  surge  level,  wave  height,  and  duration,  the  Kriebel  model  is 
recommended  as  being  the  most  suitable  for  use  by  FEMA. 
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PART  VI:  A  SYSTEMATIC  APPROACH  TO  V-ZONE  DETERMINATION 


108.  This  section  describes  a  methodology  for  inclusion  of  the  erosio:. 
process  in  V-zone  determination  which  combines  use  of  the  Kriebel  and 
Vellinga  models  with  field  measurements,  interpretation  of  historic  data,  ar.d 
engineering  judgments.  If  correctly  implemented,  the  described  procedure 
should  result  in  rational  and  consistent  V-zone  widths.  Areas  for  which  the 
procedure  is  most  applicable  are  also  identified. 

109.  It  is  worthwhile  to  note  that  the  recommendations  to  be  made 
nere ,  though  independently  determined,  are  similar  to  procedures  currently 
implemented  by  the  State  of  Florida  in  establishing  a  Coastal  Construction 
Control  Line  which  is  based  on  the  storm  surge,  waves,  and  erosion  caused  by 
a  100-year  storm  (Chiu  and  Dean  1984).  The  use  of  a  highly  simplified  ver¬ 
sion  of  the  Kriebel  model  is  augmented  by  an  extensive  program  of  field  data 
collection  which  includes  beach  and  nearshore  surveys  around  the  state.  Data 
collected  as  part  of  this  program  are  used  to  further  verify  and  improve 
erosion  predictions.  A  number  of  major  storms  have  already  been  documented 
and  are  reported  by  Chiu  (1977)  and  Balsiliie  (1985c,  1985d). 

Use  of  the  Model 

110.  We  recommend  a  straightforward  application  of  the  Kriebel  model, 
as  specified  in  the  appropriate  references.  Refinements  discussed  by  Kriebel 
(1986),  though  not  tested  here,  should  be  included.  The  Vellinga  model  will 
probably  not  be  as  widely  applicable  in  its  present  form  as  the  Kriebel  model 
but  may  be  useful  in  some  areas.  Use  of  the  MSBWT  model  is  not  yet  recom¬ 
mended  because  of  its  lack  of  independent  verification.  The  easily  computed 
erosion  estimates  given  by  Equation  13  may,  however,  be  useful  as  an 
independent  check  on  above-MSL  erosion  estimates.  Further  evaluation  of  the 
MSBWT  model  is  suggested  particularly  with  respect  to  its  simulation  of  such, 
features  as  bar  formation  and  overwash. 

111.  Actual  beach  survey  data  are  used  to  define  the  slopes  and  berm, 
features.  The  equilibrium  slope  equation  of  the  offshore,  based  preferably 
on  actual  surveys  or  grain  size  (see  Balsiliie  1982),  is  determined.  Then 
the  Kriebel  model  is  run  using  the  100-year  surge  level  (including  the  tide) 
and  duration.  FEMA  currently  computes  both  the  surge  level  and  hydrograpi. 


for  each  synthetic  hurricane  generated  for  the  determination  of  the  100-year 
surge  level;  however,  the  hydrographs  are  not  saved.  Since  for  northern 
states  the  surge  and  duration  associated  with  a  northeaster  may  be  higher  and 
longer,  respectively,  these  quantities  must  also  be  determined. 

112.  Once  the  model  has  been  used  to  compute  an  estimate  of  the  volume 
loss  above  the  surge  level,  this  estimate  is  multiplied  by  the  variability 
factor  (e.g.,  2.0  from  Part  I,  or  another  representative  value).  If  inade¬ 
quate  dunes  exist  to  survive  this  loss,  the  back-dune  area  is  allowed  to 

f lood . 

113.  Because  the  flooding  will  be  caused  by  localized  dune  breaching 
and  overwash  (unless  the  surge  level  erodes  or  overtops  the  dune  crest),  the 
depth  of  flooding  is  difficult  to  predict.  If  the  surge  is  of  short  dura¬ 
tion,  the  flooding  will  occur  as  weir  flow,  and  there  will  be  a  drop  in  water 
level  across  the  eroded  dune.  However,  once  the  dune  overwashes  inland,  the 
maximum  surge  level  should  be  assumed.  Since  the  overwash  process  is  poorly 
documented  and  unquantified,  a  conservative  and  recommended  approach  is  to 
use  the  maximum  flood  level  if  the  dune  is  flooded. 

114.  Finally,  the  results  should  be  evaluated  (i.e.,  given  a  final 
engineering  review)  for  reasonability  relative  to  historic  data  and  results 
for  adjacent  or  similar  areas. 

Historic  Data  and  Long-Term  Erosion 


115.  Considering  our  basic  lack  of  understanding  of  the  beach  erosion 
process,  the  best  and  most  direct  estimate  of  erosion  during  a  possible 
future  storm  is  made  through  use  of  measurements  collected  during  similar 
storms.  This  procedure  is  seldom  practical.  Either  the  data  are  of  insuf¬ 
ficient  quality  or,  more  frequently,  the  design  storm  has  never  occurred  or 
been  measured.  The  large  variability  discussed  in  Part  I  further  complicates 
the  situation.  However,  historic  shoreline  position  data  (surveys,  air 
photos,  high-water  marks,  etc.)  should  always  be  used  as  a  check  of  numerical 
predictions.  A  good  example  of  the  potential  use  of  historic  data  is  in  the 
determination  of  overwash  zones.  In  many  areas,  poststorm  aerial  photographs 
have  documented  overwash  penetration  limits.  Since  many  storms  are  less 
severe  than  a  100-year  event,  these  photographs  provide  a  d i rect  measure  of 
at  least  minimum  overwash  penetration  distances.  Overwash  features  are 


usually  detectable  long  after  they  are  formed.  V-zone  limits  should  ignore 
the  existing  dune  (use  the  maximum  surge  level)  in  areas  of  historic  overwash 
unless  the  shoreline  has  been  significantly  altered  through  the  construction 
of  a  substantial  duneline  or  high  seawall. 

116.  Although  long-term  erosion  is  not  of  primary  concern  during  a 
100-year  event,  its  importance  to  storm-related  losses  cannot  be  discounted. 
Quite  simply,  areas  prone  t '  long-term  erosion  will  sustain  greater  damage 
during  the  100-year  event  and  during  every  lesser  event  that  precedes  it. 
Moreover,  the  protection  offered  by  a  substantial  duneline  may  not  survive 
until  the  100-year  storm. 

117.  Long-term  erosion  rates  are  usually  based  on  measurements  of  the 
shoreline  position  inferred  from  aerial  photographs  (or  occasionally  based  on 
actual  surveys)  taken  10  to  AO  years  apart.  Usually,  only  a  few  measurements 
through  time  are  available,  so  the  rates  do  not  reflect  rapid  seasonal  or 
storm  changes.  They  do  provide  an  indication  of  how  one  section  of  shoreline 
changes  relative  to  another.  Unfortunately,  inclusion  of  a  long-term  erosion 
rate  requires  specification  of  a  time  span.  Predicted  shoreline  positions 
might  be  computed,  for  example,  based  on  a  30-year  mortgage  life  or  a  70-year 
structure  life.  Whereas  arguments  can  be  made  for  both,  use  of  a  particular 
time  span  may  be  unrealistic  in  determining  V-zones.  Yet,  some  provision  for 
recegr.  izing  high  erosion  areas  must  be  included. 

118.  Although  this  study  does  not  specifically  address  this  problem, 

it  may  be  possible  to  indirectly  include  some  measure  of  long-term  erosion  by 
adjusting  the  value  of  the  variability  factor  for  known  high  erosion  arras. 
Because  the  variability  factor  is  based  on  longshore  variability  during 
storms,  it  is  reasonable  to  assume  that  it  also  reflects  long-term  longshore 
variations  since  they  result  from  the  integrated  effect  of  many  storms. 
Therefore,  use  of  a  variability  factor  based  on  the  extreme  value  (instead  of 
the  75  percent)  could  be  used  in  high  erosion  areas  which  could  be  classified 
into  generaL  erosion  categories  of  low,  medium,  and  high  rates  based  on 
long-term  measurements.  Dunes  in  high  erosion  zones  would  have  to  be  more 
substantial  than  those  in  stable  or  low  eroding  areas.  Further  research  on 
this  or  an  alternative  method  and  a  thorough  evaluation  using  field  data  are 


Coastal  Areas  Suitable  for  Model  Application 

li9.  Three  typ^s  of  coastlines  have  been  classified  according  to 
suitability  for  model  application  along  the  5,600  km  of  the  Gulf  of  Mexico 
and  Atlantic  Ocean  coasts.  These  classifications  are  listed  in  Table  6. 


Table  6 

Classification  of  Coast  for  Application  of  Dune  Erosion  Model 


Type 

Description 

Examples 

Percent  Gulf  ana 
Atlantic  Coast* 

I.  Mostly 
suitable 

Sand  beaches  with¬ 
out  structures  or 
high  erosion;  with 
or  without  dune- 
1  Ine. 

South  Texas;  Florida 
panhandle;Outer  Banks, 
N.C. 

54 

II.  Marginally 
suitable 

Limited  extent  of 
sand  beaches, 
high  erosion 
and/or  structured. 

Portions  of  New  Jersey 
coast,  southern  Maine 

24 

III.  Not 

required 

Typically  rocky, 
cllffed,  or  other 
nonsand  beaches. 

Eastern  Long  Island, 
N.Y.;  Ten  Thousand 
Islands,  Fla. 

22 

*  Percentages  are  approximate. 


Approximately  22  percent  of  this  coast  is  nonsandy  (mud  or  rock)  where  the 
model  is  not  required  and  existing  FEMA  procedures  should  be  used.  These 
regions  are  classified  as  Type  III.  Of  the  remaining  78  percent,  approx¬ 
imately  69  percent  consists  of  sandy  beaches  where  use  of  the  model  is 
classified  as  Type  I,  "mostly  suitable."  The  remaining  coast  is  classified 
as  Type  II,  "marginally  suitable,"  either  because  of  high  annual  erosion 
rates  (greater  than  3  m/year)  or  because  of  a  high  concentration  of  coastal 
structures.  Figure  23  and  Table  7  identify  the  areas  which  are  classified  as 
either  marginally  unsuitable  or  not  required. 

120.  Within  the  Type  I  coast,  the  limitations  of  the  model  restrict 
its  use.  Obviously,  the  model  should  only  be  applied,  and  is  only  required, 
in  areas  where  the  dune  crest  is  higher  than  the  surge  level.  The  Kriebel 
model  will  function  on  low  dunes  but  only  up  to  the  point  in  the  surge  hydro¬ 
graph  when  the  dune  is  overtopped,  a  feature  which  may  be  useful.  Low  areas 
without  dunes  should  be  treated  as  flooded  profiles. 

121.  Additional  care  must  be  exercised  in  using  the  models  in  a  Type 
II  zone.  If  the  basic  assumptions  of  the  model  are  satisfied,  the  results 
should  be  meaningful.  Potential  problem  areas  are  those  where  strong 
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r r  »  ■ 


Figure  23.  Applicability  of  the  Vellinga  and  Kriebel  dune 
erosion  models  along  the  east  and  gulf  coasts  (areas  not 
marked  are  Type  I,  mostly  suitable) 


longshore  gradients  in  sediment  transport  are  known  to  exist,  such  as  near 
tidal  inlets  and  coastal  headlands. 


Role  of  E r. g  i  n e e  r  ij'.jj? J < i d gme r. t 


122.  Although  it  is  desirable  to  have  a  tool  which  would  pro  lure 
accurate  and  consistent  predictions  of  dune  erosion  and  result  ing  V-  .tor« 
widths,  no  such  model  exists  nor  is  expected  to  exist  in  the  f<  re.iee.ihl  > 


v-v-v-Vo  ->-;  v 
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Table  7 

location  and  Description  of  Mostly  Unsuitable  and  Marginally 
Suitable  Coastal  Areas 


Area  Type 


location 


Description 


1 

II 

Galveston  Island,  Tex. 

Structured 

2 

II 

Cameron  and  Vermilion  Parishes,  La. 

Eroslonal 

3 

II 

Terrebonne  Parish,  La. 

Eroslonal 

4 

II 

Terrebonne- Jefferson  Parishes,  La. 

Eroslonal 

5 

III 

Mississippi  river  delta  front,  La. 

Marsh/mud  coast 

6 

II 

Chandeleur  Islands,  La. 

Eroslonal 

7 

III 

Lighthouse  Point  to  Anclote  Key,  Fla. 

Marsh  coastline 

8 

11 

Pinellas  County  beaches,  Fla. 

Structured 

9 

II 

Anna  Marla  Island,  Fla. 

Highly  eroslonal 

10 

III 

Cape  Romano  to  Flamingo,  Fla. 

Mangrove  coast 

11 

III 

Florida  Keys  to  Coral  Gables,  Fla. 

Limestone  Islands 

12 

II 

Miami,  Fla. 

Structured 

13 

II 

North  Palm  Beach  to  Lake  Worth,  Fla. 

Structured 

14 

II 

Cocoa  Beach  to  Indian  Harbor  Beach,  Fla. 

Structured 

15 

II 

Cape  Canaveral,  Fla. 

Eroslonal 

16 

II 

Duval  County,  Fla. 

Eroslonal 

17 

II 

Hunttcg  and  Pritchards  Islands,  S.C. 

Eroslonal 

18 

II 

Capers  and  Morris  Islands,  S.C. 

Eroslonal 

19 

II 

Virginia  Beach,  Va. 

Structured 

20 

II 

Virginia  Barrier  Islands,  Va. 

Eroslonal 

21 

II 

Ocean  City,  Md.,  to  Bethany  Beach,  De. 

Structured 

22 

II 

Cape  May  to  Sea  Isle  City,  N.J. 

Structured 

23 

II 

Ocean  City  to  Seaside  Heights,  N.J. 

Structured 

24 

II 

Point  Pleasant  to  Sandy  Hook,  N.J. 

Structured 

25 

II 

Rockaway  Beach  to  Oak  Beach,  N.Y. 

Structured 

26 

II 

Southampton  to  Montauk  Point,  N.Y. 

Cl  If  fed  coast 

27 

II 

Rhode  Island,  southern  Mass. 

Rocky  or  cllffed 

28 

II 

Eastham  to  Race  Pt.,  Cape  Cod,  Mass. 

Sand/tl 1 1  cliff 

29 

II 

Plymouth,  Mass.,  to  New  Castle,  N.H. 

Rocky  or  cllffed 

30 

III 

North  of  N.H. /Mass  state  line 

Rocky  coast 

future.  Our  inability  to  understand  and  correctly  model  the  underlying 
processes,  combined  with  the  natural  variation  and  the  uncertainty  in 
predicting  future  storms,  precludes  the  development  of  a  model  based  on  fir 
principles.  Consequently,  we  must  use  simpler  models  for  guidance  and 
combine  that  guidance  with  other  sources  of  data,  as  described  above,  and 
experienced  engineering  judgment. 

12'.  Similar  guidance  is  a  cornerstone  of  other  FEMA  procedures  Us 
Of  historical  data  and  engineering  judgment  is  an  often  repeated  recommer, 
da: in  in  the  "Guidelines  for  Identifying  Coastal  High  Hazard  Zones"  (US  Ar 
Kr.p  ;  r.eer  I.’is'ri'-t  ,  Galveston  I97S ) .  With  reference  to  sand  dunes,  this 

r  e  ;  S  r  t  state:;, 

In  many  i  nut  ,ir.--es  ,  r;ul  .tantial  sand  dunes  with  heights  of  1  5 
f *  at  ve  sea  level  and  base  widths  of  several  hundred  feet 
have  been  swept  away  by  storm  surge  and  a<  company  i  ng  wave:;. 

Unless  histori'-al  data  indicate  t  hat  sand  dunes  in  the  area 
have  repeatedly  withstood  wave  attack  during  storms,  they 
sh-.uld  n.,{  go  considered  as  effect  ive  surge  and  wave  hart  iers. 


Although  presently  available  models  can  objectively  help  to  reduce  the 
uncertainty,  this  advice  remains  true  today.  When  the  engineer  is  faced  with 
a  marginal  duneline,  it  should  not  be  considered  an  effective  wave  barrier. 


PART  VII:  SUMMARY  AND  CONCLUSIONS 


124.  This  report  has  examined  the  merits  of  both  theoretically  based 
and  parameterized  models  for  possible  use  in  estimating  the  impact  of  dune 
erosion  in  mapping  high  hazard  V-zones.  Among  these,  only  the  empirica  :  iy 
based  models  were  shown  to  produce  reasonable  results  with  a  minimum  of 
effort  and  required  input  data.  Of  the  models  evaluated  in  this  feasibility 
study,  the  Kriebel  model  was  judged  superior  to  the  Vellinga  model,  not 
because  it  produced  more  accurate  results  but  because  it  (a)  is  applicable  to 
a  wider  range  of  beach  configurations;  (b)  is  less  demanding  in  terms  of  the 
accuracy  requirements  of  the  input  data:  and  (c)  accounts  for  storms  of 
varying  duration.  This  last  capability  is  particularly  important  because 
Atlantic  coast  beaches  are  affected  by  both  short  duration  hurricanes  and 
long  duration  northeasters. 

125.  Currently,  the  Kriebel  model  does  not  fulfill  all  the 

requi rement s  of  an  ideal  erosion  model.  Its  limitations  include  an  inability 
to  model  offshore  bars,  a  dependence  on  an  overly  simplistic  profile  shape, 
and  an  insensitivity  to  varying  wave  conditions.  Nevertheless,  if  properly 
implemented,  the  Kriebel  model  is  expected  to  produce  reasonably  accurate 
dune  erosion  estimates,  and  it  is  believed  that  the  capability  of  the  model 
may  be  improved  in  the  future. 

126.  The  foregoing  conclusion  is  based  on  examination  of  available 
models,  including  the  comparison  given  in  Part  IV.  This  comparison  augmented 
the  limited  verification  of  both  the  Kriebel  and  Vellinga  models  which  had 
been  dene  to  date.  The  Kriebel  model  was  originally  developed  as  a  Master's 
thesis  ar.d  was  evaluated  with  a  single  Florida  profile  surveyed  after  the 

117  5  Hurricane  El  >i se .  Kriebel  (1936)  provided  additional  verification  us i r g 
20  Florida  profiles  from  Hurricane  Eloise.  Considerable  data  now  exist  which 
would  perr.it  a  more  extensive  evaluation.  For  example,  the  State  of  Florida 
has  hundreds  of  profiles  surveyed  after  storms  (Balsillie  1985c,  d;  Chiu 
1  9  7  7  )  . 

127  .  Time  constraints  of  this  study  permitted  only  a  limited  a:;.  ;r.t  .  ; 
test  ing  t  <■  be  conducted .  Consequent  I  y  ,  if  the  recommended  Kriebel  n.  del  ; 
adopted,  the  •  .<per i enoe  that  FEMA  gains  through  use  of  the  model  on  a  wide; 
variety  of  coastal  areas  should  be  applied  to  improving  the  model  ,,r,  !  :  t  s 

irr.pl  err.®  nt.it  ion. 


i 


128.  The  reviewed  MSBWT  model  has  a  number  of  potentially  attractive 
features,  including  use  of  a  probabilistic  approach  to  erosion  quantities, 
provision  for  overwash,  and  a  duration  dependency  based  on  the  rise  time  of 
the  peak  surge.  Though  it  is  too  complex  and  undeveloped  for  the  present 
needs  of  FEMA ,  this  model  may  be  useful  in  the  future. 

129.  Because  none  of  the  models  can  adequately  account  for  the  norma 
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APPENDIX  A:  DISCUSSION  AND  ANNOTATED 


BIBLIOGRAPHY  OF  OVERWASH  LITERATURE 


I  nt  roduc  t  i  on 


1.  Overwash  is  defined  by  the  American  Geological  Institute  (Bates  and 


Jackson  1980)  as: 


A  mass  of  water  representing  the  part  of  the  uprush  that  runs 
over  the  berm  crest  (or  other  structure)  and  that  does  not 
flow  directly  back  to  the  sea  or  lake. 


The  flow  of  water  in  restricted  areas  over  low  parts  of 
barriers  or  spits,  especially  during  high  tides  or  storms. 


Although  this  definition  applies  to  any  low  shoreline,  including  reef  islands 


>r  rocky  spits,  in  this  study  only  overwash  along  sandy  barrier  beaches, 


especially  if  a  frontal  duneline  is  present,  is  discussed. 


A  washover  is  the  geomorphic  feature  produced  by  an  overwash  event. 


The  passage  through  the  duneline  is  called  the  throat  (Figure  Al).  The  fan 


is  the  portion  of  the  washover  deposited  on  the  island  or  beach  interior 
where  the  landward  flow  is  unrestricted  by  topography  (Figure  Al).  Effects 


of  coastal  storms  on  beaches  are  progressive.  Initial  beach  erosion  and  dune 


scarping  will  be  followed  by  slight  overtopping  of  the  lower  portions  of  the 


duneline  as  the  water  level  and  wave  runup  increase.  With  continued  water 


level  rise,  significant  dune  overtopping  will  occur,  with  associated  fan 


deposition  (Figure  A2 ) .  If  the  surge  level  rises  enough,  complete  over¬ 


topping  of  the  duneline  will  occur,  possibly  resulting  in  total  erosion  of 


the  frontal  dune.  Significant  storm  duration  without  an  increase  in  surge 


.evel  may  produce  a  similar  progression  of  effects,  with  dune  erosion  causing 


.'erwasn  to 


occur  at  lower  water  levels.  If  landward  flow  of  water  is 


t  rated  in  o.ne  area,  deep  erosion  of  an  overwash  throat  may  occur, 


!;  i:.,-  in  formation  of  a  tidal  inlet,  although  most  tidal  inlets  formed 


of  elevated  bav  water  (Hite  1 0 2 -I 


Hurricane  of  21  September  1938  upon  the  north  Atlantic  coast  from  Long 
-sland.  New  York,  to  Cape  Cod,  Massachusetts,  were  described  by  several 
researchers  (Wilby  et  al.  1939  and  Howard  1939).  The  storm  made  landfall  in 
southwest  Rhode  Island,  with  an  estimated  tide  of  4.9  m  at  Point  Judith 
:Wiiby  et  al.  1939).  Although  this  storm  was  of  short  duration,  with  the 
rise  and  fall  of  the  storm  force  ending  in  5  to  6  hr,  total  destruction  of 
coastal  con.mur.it  ies  occurred.  The  significance  of  overwash  in  widening  the 
islands  and  importance  of  storm- induced  inlet  formation  was  noted  by  these 
authors.  Wilby  et  al.  ; 1939)  concluded  that  barrier  dunes  provide  some 
protection  tor  coastal  buildings  during  storms,  and  wider  beaches  may  afford 
protect  ion  to  landward  areas. 

4-  Mar.y  of  the  earlier,  qualitative  studies  of  overwash  were  completed 
alo.ig  the  gulf  coast  (Morgan  1959,  Fisk  1959,  Scott,  Hoover,  and  McGowan 
1969,  and  Hayes  1967).  In  this  microtidal  setting  (tidal  range  less  than 
“  m'>  barriers  are  typically  long  and  linear,  with  few  permanent  tidal 
inlets.  Deposition  of  continuous  overwash  aprons  (coalescing  washover  fans) 
occurs  on  the  landward  sides  of  barriers  during  storm  passage.  More  recent 
work  has  revealed  that  complete  overwash  of  low  profile  Gulf  barriers  may 
occur  during  hurricane  passage  (Otvos  1979  and  Nummodal  1982),  resulting  in 
pos.  storm  submergence.  Years  later,  emergent**  and  redevelopment  m.iy  neons  in 
.more  landward  positions,  as  documented  by  Otvos  (19, '9)  for  portions  of  the 
Char.de  lour  Island  chain  of  Louisiana. 

5.  Cne  o»  the  tirst  detailed  analyses  of  a  ..ingle  w  a  shove;  fan  was 
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sediment  contribution  from  various  processes.  A  summary  of  the  contribution 
of  landward  sediment  transfer  processes  at  several  Atlantic  coast  localities 
is  contained  in  Table  Al.  Eighty-five  percent  of  the  North  Carolina  coast 
has  been  subjected  to  overwash  between  1938  and  1974  (Boc  and  Langfelder, 
1977).  These  data  indicate  that  the  magnitude  of  overwash  varies 
significantly  with  locality. 


Table  Al 


Locality 

Overwash  Wind 

X  X 

Delta 

X 

Data  Source 

Outer  Banks, 

North  Carolina 

13 

14 

73 

Pierce  (1969) 

Assateague  Island 
Maryland 

,  12 

6 

62 

Bartberger  (1976) 

Malpeque  barrier 
system,  Canada 

48* 

52 

Armon  (1979) 

Rhode  Island 

43* 

57 

Fisher  and  Simpson 

(1979) 

*  Cont>1  ned  effects  of  wind  and  overwash. 

7.  A  detailed  study  of  overwash  texture,  composition,  and  structure 
was  completed  by  Schwartz  (1975).  Volume  estimates  from  individual  storms 
were  also  compiled.  During  a  northeaster  ( extrat ropical  storm)  in  February 
1973,  a  minimum  estimate  of  302,000  of  washover  occurred  between  Caffeys 
Inlet  and  Buxton,  North  Carolina  for  an  average  washover  volume  of  2.7  m^/m 

t  beach f  runt  for  the  1 1 3-kra  stretch  of  coastline.  Assuming  that  only  20  to 
2  5  percent  of  t  he  barrier  experienced  overwash,  the  volume  of  overwash  per 
unit  of  overtopped  beach  was  10.78  to  13.5  m Vm  of  beach  (Schwartz  1975). 

8.  Actual  field  measurement  of  overwash  processes  was  accomplished  ty 
richer  ,  I. eat  her  man ,  and  Ferry,  (1974)  at  Assateague  Island.  Overwash 

ve  ;  .  ..•  ;  es  we  ;  <•  measured  during  st  urn  events,  and  ac/Mir  at  e  measurements  were 

"lie  »  am  ,n*  •.  tea-  *  er  >n  ,  dice  s  ai;  ;  ng  ,  and  uvet  wash  -de;  s  ;  •  ;  n 
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a  rr.e.tr,  f  1  .  S  it.  ,  e  F.  v  dept  hr.  averaged  0 . 2  m . 
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’.e  ;  p,’  i  .  r  he  rr.  at  imputant  parameter  in  let  on:.;  r.  .  ng  t  tie  ii.a  gi,  :  t  ude  ,  1  an 
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overwash.  During  a  26-month  study  along  Assateague  Island,  seven  overwash 
events  were  recorded.  Although  backshore  deposition  occurred  during 
overwash,  a  slight  decrease  in  volume  occurred  because  of  poststorm  eoiian 
transport  back  to  the  beach  (Leatherman  1976).  Net  seaward  transport  of 
overwash  deposits  by  wind  was  also  documented  along  the  Tabusintac  barrier 
system  in  New  Brunswick,  Canada  (Rosen  1979). 

10.  The  influence  of  human  activity  upon  natural  beach  processes  has 
been  a  subject  of  study  for  many  years.  One  major  debate  involves  the 
effects  of  artificial  dunes  upon  storm-induced  beach  erosion.  Some 
researchers  believe  that  increased  turbulence  resulting  from  wave  backwash 
frcm  the  dune  front  results  in  increased  beach  erosion.  These  researchers 
also  believe  that  an  artificial  duneline  will  prevent  island-building 
oV'Twas!: .  Other  researchers  believe  that  erosion  of  artificial  dunes  during 
storms  nourishes  eroding  beaches,  protecting  against  more  severe  beach 
or  •;  i  on  ( Leatherman  1981).  Other  human  activities  influencing  natural  storm 
;r  -oases  include  the  formation  of  overwash  channels  at  locations  of  teach 
a  s  cuts  through  the  duneline  which  occurred  along  the  eastern  port  ion  cf 
Os  :■■■••  - 1  .  n  Island,  Texas,  during  Hurricane  Alicia  (Dupre  1985  ).  Although 
s. :  • :  r.x  ar»  often  used  to  stabilize  beaches,  during  severe  storms  they  rr.av 
:  r  '  r «  2  i  i;ted  erosion  and  overwash,  as  shown  in  Figures  A3a  and  b  :  "S 
Arr.v  engineer  District,  New  Orleans  1965). 
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sediment  that  had  been  carried  across  a  barrier  island  and  dumped  in  the 
adjoining  bay  by  catastrophic  storms,  such  as  hurricanes.  One  such  fan, 
measuring  4-1/4  by  4  miles  (6.9  by  6.4  km),  occurs  at  the  north  end  of  St. 
Joseph  Island,  on  the  central  Texas  coast.  It  grades  northeastward  into  a 
sutaerial  tidal  delta.  The  surface  facies  of  these  two  features  and  the 
adjoining  barrier  island  nucleus,  as  well  as  the  geometry  of  the  component 
sediment  bodies,  are  described  in  detail.  The  data  upon  which  description 
and  interpretations  are  based  were  obtained  by  digging  trenches  and  pits, 
to  6  ft  (2.4  m )  deep. 

The  washover  fan  is  fed  via  a  low-level  washover  channel,  Vincent 
Slough,  which  is  cut  through  the  barrier  island  nucleus.  The  fan  is  accre 
ir.g  laterally  into  the  adjoining  bay,  Aransas  Bay,  through  the  deposition 
she.!  debris  and  sand  that  are  derived  from  the  inner  neritic  zone  of  the 
':.t  ;  rent  a  1  shelf .  Periodic  hurricane  surges,  both  flood  and  ebb,  deposit 
large  volumes  of  sediment  in  Aransas  Bay.  The  washover  fan  thus  built  up 
1 “  to  1 0  in.  (63.5  to  76.2  cm)  thick  at  its  periphery  and  thickens  to  more 
t  hur.  r  J  in.  I  127  cm)  near  its  apex.  It  is  composed  predominantly  of  an 

r  i  -ate  set  of  sheet  sands,  each  3  to  15  in.  !  7 . 6  to  3  8.1  cm)  thick.  Eu 
w  !  v r  sand  rests  r.  a  sharp  scour  base  and  consists  of  a  shell  placer  a 
•  r  i  e  '..at  grades  upward  into  horizontally  laminated,  shell-free,  clean 
d  . 

The  washover  fan  is  a  two-part  feature.  One  part  is  still  actively 
:  :  .  :  a  ; : nr ,  albeit  very  slowly,  and  mainly  consists  of  sand  of  the  w a  shove 
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r. ,  jr.j  Simpson ,  E.  J.  1979.  Washover  and  Tidal  Sedimentation 
-  .is  Fr.v  i  r  omental  Factors  in  Development  of  a  Transgressive  Barrier 

:  inline,"  in  L.ea  t  he  rni.i  r.  ,  S.  P.,  ed .  ,  Barrier  Islands:  From  the  Gulf  of  St 

:  •  ■  the  f  of  Mexico ,  Academic  Press,  New  York,  N.Y.,  pp  127-148 

Washover  fans  and  tidal  deltas  are  significant  reservoirs  of  sediment 

:  : er  island  systems.  These  backbarrier  and  lagoonal  deposits  also 

:  e;  :  eser.t  important  sediment  sinks  in  the  overall  littoral  sediment  budget 

a  s:..  reline.  Along  the  40-km-long  Rhode  Island  south  shore  barrier  beaches 

ex:  e  :  ve  wash,  over  fans,  resulting  from  hurricanes,  and  several  inlet 

its  ate  prevalent  .  A  long-term  quantitative  analysis  of  the  subtidal 

.;  rat  idal  sedimentation  of  these  units  from  1939-  1975  was  conducted 

;  s:  r.  g  t  o  rammot  r  i  s  techniques.  Backbarrier  sedimentation  features  were 

i  dent  i  f  ted  in  th.e  field  to  develop  "ground-truth"  keys  to  identify  it. 

it  .res  .  n  t  he  aerial  photography.  Sediment  units  were  measured  on  each 

.  f  aer  ;  a  1  ;  -t  •.•graphs  by  point-counting  of  grids.  These  photographic 

:•  •’  asurerr.er.t  s  were  scaled  by  using  ground  survey  data. 

T.ring  the  3  0 -year  study  period,  the  total  areal  sedimentat  ion  change 

su;  rati  d  a  1  washover  deposits  was  +522,790  m^  .  The  subtidal  washover 

detnsits  am.,  anted  to  +267,950  ,  for  a  total  washover  sedimentation 

o 

accumulation  of  790 , 740  m-  .  For  the  same  time  period,  the  total  area  chin.- 
of  supratidai  flood  tidal  deltas  at  the  inlets  was  +188,240  m~  ,  while  r 
the  subtidal  tidal  deposits  was  +862,320  for  a  total  tidal  delta 
deposition  of  +1,050,560  m^  .  This  analysis  indicates  that  during  this 
36-year  period,  tidal  delta  sedimentation  is  1-1/3  times  more  effe  • 

washover  in  ‘.he  landward  transportation,  deposition,  and  storage  f  • 

This  barrier  coast,  in  addition  to  transgressing,  is  ere- 1  i  f.g  •  •  • 

of  0.7  m/yr.  An  inverse  relationship  exists  between  overwash  •••■■ 

barrier  island  width,  which  can  be  further  related  to  be  a 
Measured  beach  erosion  ever  the  36-year  period  is  dire-r  lv  : . 
rate  of  overwasn  occurrences  at  different  points  ai  r..- 


Fisher,  J.  S.  1980.  "Field  and  Laboratory  Sr  .  ry 
Dynamics,"  Final  Report,  US  Army  Resear  -h  ■ 

Field  studies  at  Asset  eague  I  s  1  a  :r  ’  ■  ’  ' 

overwash.  Repeated  surveys  at  sn'- :  i  .  w  = 
generated  overwash  transports  s  -  •  * 

for  storms  characterized  by  ;  • 

this  material  seaward.  ;  .< . 
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Godfrey,  P.  J.,  Leatherman,  S.  P.,  and  Zaremba,  R.  1979.  "A  Geobotanical 
Approach  to  Classification  of  Barrier  Beach  Systems,"  in  Leatherman,  S.  P., 
ed . ,  Barrier  Islands:  From  the  Gulf  of  St.  Lawrence  to  the  Gulf  of  Mexico, 
Academic  Press,  New  York,  N.Y. ,  pp  99-126. 

Barrier  beaches  from  Cape  Cod  to  Cape  Lookout  have  been  studied  using 
geological  and  ecological  techniques.  The  resulting  data  have  beer,  inte¬ 
grated  to  demonstrate  the  interrelationships  among  plants,  processes,  and 
barrier  morphology.  From  this  analysis  a  general  pattern  for  east  coast 
barriers  undergoing  recession  emerges.  The  regional  variation  in  vegetation 
and  its  response  to  overwash  have  been  found  to  be  an  important  criterion  for 
classification.  The  vegetation  of  Northern  barrier  beaches  is  dominated  by 
American  beach  grass  (Ammophila  breviligulata)  in  the  dune  strand  community 
and  by  the  decumbent  form  of  salt  meadow  cordgrass  (Spartina  patens)  in  the 
high  marsh.  Salt  meadow  cordgrass  is  killed  by  overwash  burial  and  is 
replaced  by  dune  vegetation,  originating  from  seeds  and  plant  fragments  found 
in  the  drift  lines.  The  stratigraphy  of  a  transgressive  barrier  in  the 
Northeast  shows  a  sharp  demarcation  between  the  salt  marsh  and  overlying 
washover /dune  sands. 

In  the  Southeast,  sea  oats  (Uniola  paniculata )  dominate  the  dune 
grasslands,  while  the  upright  variety  of  Spartina  patens  is  ubiquitous.  The 
dunes  along  the  Outer  Banks  develop  initially  as  scattered  clumps  resulting 
from  lack  of  well  developed  drift  lines,  irregular  seed  dispersal,  and  the 
clumped  growth  deposition  on  barrier  flats  and  marshes  in  the  absence  of  a 
continuous  barrier  duneline.  The  upright  form  of  Spartina  patens  has  the 
ability  to  grow  through  this  overwash  sediment  and  reestablish  itself  on  the 
fan  surface.  An  analysis  of  sedimentary  sequences  shows  a  deposit  of  clean 
overwash  layers  alternating  with  organic  layers.  These  biogeological  studies 
have  shown  the  importance  of  regional  variation  in  the  vegetation  in  deter¬ 
mining  barrier  beach  topography.  (Authors) 

Hayes,  M.  0.  1967.  "Hurricanes  as  Geological  Agents:  Case  Studies  of 

Hurricanes  Carla,  1961,  and  Cindy,  1963,"  Texas  Bureau  of  Economic  Geology, 
Report  of  Investigations  No.  61. 

Tropical  storms,  which  cross  the  Texas  coastline  with  a  frequency  of 
0.67  storms  per  year,  play  a  major  role  in  nearshore  sedimentation  on  the 
south  Texas  coast.  Greatest  geological  effects  of  these  storms  are  produced 
by  wind-driven  waves  and  by  storm  surges. 

The  comparison  of  a  part  of  the  nearshore  environmental  complex  of  a 
segment  of  the  south  Texas  coast  before  and  after  Hurricane  Carla,  1961, 
shows  the  effects  of  the  storm.  The  storm  removed  a  belt  of  foredunes  20  to 
50  yards  (15  m)  wide  from  the  seaward  side  of  Padre  Island  and  left  the 
foredune  ridge  with  wave-cut  cliffs  up  to  10  ft  (3  m)  high.  The  formation  of 
a  broad,  flat  hurricane  beach  drastically  altered  the  beach  profile.  The 
landward  side  of  the  barrier  island  (wind-tidal  flats)  received  much  washover 
material  containing  surf  zone  and  beach  mollusks.  The  storm  also  submerged 
high-level  mud  flats  along  the  landward  side  of  Laguna  Madre  and  covered  them 
with  a  fresh  layer  of  mud.  A  much  milder  storm  (Cindy)  passed  through  the 
area  in  September  1963,  and  a  small  swash  bar  was  deposited  over  the  seaward 
edge  of  the  preexisting  hurricane  beach.  (Author) 


Kahn,  J.  H. ,  and  Roberts,  H.  H.  1982.  "Variations  in  Storm  Response  Along  A 
Microtidal  Transgressive  Barrier-Island  Arc,"  Sedimentary  Geology,  Vol  33, 
pp  129-146. 

Storm  response  along  the  transgressive  Chandeleur  barrier  island  arc 
southeast  of  the  Mississippi  Delta  plain  is  variable  because  of  local  dif¬ 
ferences  in  sediment  supply,  shoreline  orientation,  and  barrier  morphology. 

A  study  of  the  morphological  impact  of  Hurricane  Frederic  (1979)  affirmed 
that  tropical  storms  are  the  primary  agents  causing  erosion  and  migration  of 
this  barrier  arc. 

Frederic's  greatest  impact  was  in  the  duneless  southern  Chandeleurs 
where  sheet-flow  overwash  caused  flattening  of  the  barrier  profile,  destruc¬ 
tion  of  a  strip  of  marsh  50  to  100  m  wide,  and  shoreline  retreat  of 
approximately  30  m.  In  contrast,  overwash  in  the  northern  Chandeleurs  was 
confined  between  dunes  in  channels  established  by  previous  storms.  This 
channelized  overwash  breached  the  northern  Chandeleur  barriers  in  19  places. 
As  Frederic  passed,  return  flow  through  these  channels  transported  overwash 
sediment  back  to  the  nearshore  zone.  These  ebb  deposits  were  a  source  for 
longshore  drift  sediments  which  quickly  sealed  storm  channels,  reestablishing 
a  coherent  northern  Chandeleur  barrier  arc. 

These  storm  response  patterns  may  help  explain  long-term  changes  in 
barrier  morphology.  During  an  84-year  period  (1885-1969),  the  southern 
Chandeleurs  decreased  41  percent  in  area,  with  an  average  retreat  rate  of  9.1 
m/yr,  compared  to  a  15  percent  increase  in  area  and  an  average  shoreline 
retreat  rate  of  7.2  m/yr  for  the  northern  Chandeleurs.  (Authors) 

Kraft,  J.  C. ,  et  al.  1979.  "Processes  and  Morphologic  Evolution  of  an 
Estuarine  and  Coastal  Barrier  System,"  in  Leatherman,  S.  P. ,  ed.,  Barrier 
Islands:  From  the  Gulf  of  St.  Lawrence  to  the  Gulf  of  Mexico,  Academic  Press, 
New  York,  N.Y. ,  pp  149-184. 

Coastal  barriers  in  Delaware  are  rapidly  transgressing.  Wave-dominated 
processes  are  the  most  important  single  factor  in  determining  the  volume  of 
sediment  in  motion.  Availability  of  sediment  from  the  continental  shelf, 
erosion  of  the  barrier  itself,  and  erosion  of  headlands  determines  whether  or 
not  a  barrier  can  exist  and  evolve. 

Studies  of  the  internal  structure  of  the  barriers  indicate  that  washover 
processes  and  flood-tidal  delta  deposition  are  the  dominant  factors  in  the 
landward  transgression.  Wind-transported  sand  is  of  secondary  importance, 
derived  from  the  beach  berm  and  washover  fans  and  transported  to  coast 
parallel  dunes.  (Authors) 

Leatherman,  S.  P.  1976.  "Barrier  Island  Dynamics:  Overwash  Processes  and 
Eolian  Transport,"  Proceedings  of  the  15th  Coastal  Engineering  Conference, 
American  Society  of  Civil  Engineers,  pp  1958-1974. 

The  northern  8  km  of  shoreline  at  Assateague  Island,  Maryland,  are 
presently  being  eroded.  During  storms,  swash  surges  are  able  to  overtop  the 
most  landward  (storrr.s  berm  as  overwash,  with  deposition  occurring  on  the 
barren  flats.  Where  primary  barrier  dunes  still  exist,  sediment-charged 


surges  are  funneled  through  breaches  in  the  dune  field  for  deposition  of  the 
entrained  material  on  the  washover  fan. 

Sediment  budget  computations  show  that  there  has  been  a  small  net  loss 
of  material  at  each  washover  area,  in  spite  of  seven  discrete  overwash  events 
during  a  26-raonth  time  interval.  The  predominant  northwest  winds  effectively 
eroded  the  overwash  material,  transporting  the  majority  of  the  sand  back  to 
the  beach.  This  analysis  indicates  that  there  exists  a  balance  between 
overwash  and  eolian  processes  with  wind  transport  being  slightly  dominant. 
(Author) 


Leatherman,  S.  P.,  Allan,  T.  W. ,  and  Fisher,  J.  S.  1977.  "Overwash 
Sedimentation  Associated  with  a  Large-Scale  Northeaster,"  Marine  Geology,  Vol 
24,  pp  109-121. 

The  1  December  1974  northeaster  was  a  significant  event  in  terms  of 
sediment  transport  with  20  m^  of  sand  per  meter  of  dune  breach  being  carried 
onto  the  backdune  area  of  Assateague  Island,  Maryland,  as  overwash.  Previous 
investigators  have  reported  larger  transport  rates  for  landfall  hurricanes, 
but  this  is  the  largest  amount  recorded  for  a  winter  northeaster.  The  Ash 
Wednesday  Storm,  6-8  March  1962,  was  a  much  larger  event,  but  no  quantitative 
data  exist  for  overwash  deposition.  (Authors) 

Leatherman,  S.  P. ,  ed .  1981.  Overwash  Processes,  Benchmark  Papers  in 

Geology,  Vol  58,  Hutchinson  Ross,  Stroudsburg,  Pa. 

This  book  includes  a  collection  of  reprints  of  some  of  the  more  signi¬ 
ficant  historic  and  recent  articles  associated  with  overwash  research.  The 
book  is  divided  into  six  parts:  (1)  early  works;  (2)  coral  reefs,  islands, 
and  cays;  (3)  estuaries  and  lakes;  (4)  barrier  islands;  (5)  overwash  deposits 
in  geologic  record;  and  (6)  management  implications.  The  significance  of 
articles  within  each  section  is  discussed  by  the  editor. 


Leatherman,  S.  P.,  and  Zaremba,  R.  E.  1986.  "Dynamics  of  a  Northern  Barrier 
Beach:  Nauset  Spit,  Cape  Cod,  Massachusetts,"  Geological  Society  of  America 
Bulletin,  Vol  97,  pp  116-124. 

Quantitative  analysis  of  historical  shoreline  (ocean  and  bay)  and 
coastal  environments  (dune,  salt  marsh,  shrub,  washover,  and  sandy  beach) 
change  data  for  Nauset  Spit,  Cape  Cod,  Massachusetts,  are  described  in  this 
paper.  There  are  marked  differences  in  the  relative  role  of  the  landward 
sediment  transport  processes  along  the  Nauset  barrier  chain.  In  general,  the 
barriers  have  been  narrowing  during  the  past  century,  and  bayshore  accretion 
has  not  been  equal  to  ocean  shoreline  erosion  rates.  Inlets  are  largely 
responsible  for  bayside  sedimentation,  whereas  overwash  is  not  effective  in 
maintaining  width  until  the  barrier  slims  considerably.  While  salt  marshes 
form  on  intertidal  washover  deposits,  the  vast  expanse  of  barrier-related 
marshes  have  developed  atop  flood-tidal  deltas,  and  overwash  subsequently  has 
buried  and  killed  the  marshes.  Therefore,  massive  washovers  eventually 
become  the  sites  of  new  dunes,  often  at  the  expense  of  marshes  in  the  cyclic 
events  of  barrier  evolution.  (Authors) 


Maynard,  A.  K.,  and  Suter,  J.  R.  1983.  "Regional  Variability  of  Washover 
Deposits  on  the  South  Texas  Coast,"  Transactions,  Gulf  Coast  Association  of 
Geological  Societies,  Vol  33,  pp  339-346. 

Overwash  processes  play  an  important  role  in  determining  the 
stratigraphy  of  raicrotidal  barrier  islands.  Along  the  microtidal  coast  of 
south  Texas,  regional  variability  in  barrier  island  geometry  produces  a 
spectrum  of  washover  types.  South  Padre  Island  is  a  low-profile, 
transgressive  feature  with  a  discontinuous  to  nonexistent  foredune  ridge.  It 
displays  sheet  overwash,  coalescing  washover  terraces,  and  washover  fans  fed 
by  large  hurricane  channels.  North  Padre  and  Mustang  Islands  are  high- 
profile  barrier  islands  with  continuous  foredune  ridges,  and  thus  only  small, 
discrete  interdune  washovers  occur.  The  relict  tidal  inlets  between  Mustang 
and  Padre  Islands  are  the  sites  of  the  largest  washovers  in  the  system, 
termed  reactivated  tidal  deltas.  This  last  type,  although  relatively  rare  in 
modern  washover  deposits,  is  probably  similar  in  mode  of  origin  to  the  large 
lobate  backbarrier  features  found  on  many  high-profile  barriers  on  the  Texas 
coast.  (Authors) 


Pierce,  J.  W.  1970.  "Tidal  Inlets  and  Washover  Fans,"  Journal  of  Geology, 
Vol  78,  pp  230-234.  ~  ’ . ' 

Tidal  inlets  and  washover  fans  are  genetically  related.  The  resulting 
feature  is  dependent  upon  barrier  configuration,  depths  in  the  lagoon  adja¬ 
cent  to  the  barrier,  and  the  direction  from  which  the  storm  surge  came, 
either  from  the  sea  or  the  lagoon. 

Attack  on  barrier  islands  from  the  seaward  side  by  waves  overtopping  the 
barrier  will  result  in  washover  fans  on  wide  barriers  where  extensive  adjoin¬ 
ing  tidal  flats  are  present.  Inlets  can  be  cut  by  this  type  attack  on  narrow 
barriers  where  no  tidal  flats  are  present.  Storm  surge  from  the  lagoonal 
side,  if  channeled  along  tidal  creeks,  can  easily  cut  inlets  through  a 
barrier.  (Author) 


Rosen,  S.  P.  1979.  "Eolian  Dynamics  of  a  Barrier  Island  System,"  in 
Leatherman,  S.P.,  ed.,  Barrier  Islands:  From  the  Gulf  of  St.  Lawrence  to  the 
Gulf  of  Mexico,  Academic  Press,  New  York,  N.Y.,  pp  81-98. 

Eolian  transport  of  sand  was  monitored  by  using  directional  sand  traps 
on  the  beach,  foredune,  overwash  and  spit  environments  of  Tabusintac  barrier 
system  in  northeast  New  Brunswick.  The  total  volume  of  sand  moved  over  the 
system  was  computed  from  the  measured  transport  rates. 

A  large  volume  of  sand  (1,720  m^)  was  transported  in  the  alongshore 
directions  from  storm  overwash  deposits  into  vegetated  backdune  areas  to 
produce  vertical  accretion.  In  1977,  the  storm  overwash  deposits  accounted 
for  8  percent  of  the  subaerial  volume  of  the  system. 

Net  cross-island  transport  was  offshore  from  overwash  areas  and  the 
mid-backbeach  (1,556  m^)  in  response  to  prevalent  offshore  winds.  A 
comparable  volume  (2,100  m^)  of  sand  was  moved  onshore  from  the  beach  to  the 
foredune  base.  Approximately  220  m^  of  this  amount  crossed  accretional  dune 
crests,  but  very  little  was  transported  up  wave-eroded  dune  crests;  hence, 
only  the  former  grew  vertically.  The  remainder  was  transported  alongshore  in 
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the  direction  of  net  wave-induced  longshore  transport  and  resulted  in  the 
extension  of  the  foredune  into  overwash  and  spit  areas. 

Overwash  deposits  serve  as  a  source  of  sand,  and  these  breaches  in  the 
barrier  dune  also  act  as  corridors  for  eolian  transport.  The  redistribution 
of  storm  overwash  deposits  by  wind  is  important  in  effecting  the  vertical 
growth  of  a  barrier  island.  (Author) 


Schwartz,  R.  K.  1975.  "Nature  and  Genesis  of  Some  Storm  Washover  Deposits," 
CERC  TM-61,  US  Army  Engineer  Waterways  Experiment  Station,  Coastal 
Engineering  Research  Center,  Vicksburg,  Miss. 

Under  certain  conditions  large  amounts  of  sediment  are  moved  across  the 
beach  during  storms  and  stored  subaerially  instead  of  being  eroded  from  the 
beach  and  foredune  and  introduced,  or  simply  returned  seaward,  into  the 
littoral  drift  system.  Washover  occurrence  is  a  function  of  the  degree  of 
storm  surge  and  backshore-f oredune  relief.  The  shape  and  dimension  of  a 
washover  sand  body  are  largely  controlled  by  the  surrounding  topography.  The 
major  immediate  source  of  washover  sediment  during  storm-surge  flood  is  the 
beach  and  shoreface;  the  foredune  may  contribute  variable  amounts.  Subenvi¬ 
ronments  of  the  backbarrier  supply  sediment  during  storm-surge  ebb  flow,  and 
eolian  processes  add  variable  amounts  of  sediment  during  washover  deposition 
depending  upon  storm  wind  intensity  and  duration.  In  turn,  washover  deposits 
serve  as  a  sediment  source  to  the  poststorm  eolian  system  or  encroaching 
bodies  of  water,  e.g.,  a  transgressing  sea. 

Both  structural  and  textural  properties  are  explained  by  pulsating 
(discontinuous)  unidirectional  flow  of  sediment-charged  water  across  the 
washover  surface.  Velocities  imparted  by  the  fluid  surge  of  overwashing 
storm  waves,  combined  with  gravity  effects  due  to  a  landward-sloping  surface, 
result  in  the  development  of  flow  similar  to  sediment  gravity  flow.  Each 
flow  event  results  in  initial  scouring  of  the  subaerial  surface  followed  by 
the  development  of  planar  (horizontal)  stratification,  while  injection  into 
the  more  distal  subaqueous  setting  results  in  the  development  of  delta- 
foreset  structures.  Normal  and  inverse  textural  grading  as  well  as  textural 
coarsening  in  the  direction  of  flow  are  also  the  products  of  this  flow. 

The  observed  washover  properties  and  genetic  interpretations  describe 
washover  response  to  two  particular  storm  conditions.  Although  the  same 
basic  response  is  expected  for  storms  of  various  magnitudes,  modifications 
are  expected.  In  response  to  a  larger  storm  surge  (e.g.,  hurricane)  or  a 
different  topographical  setting,  storm  surge-ebb  may  be  an  important  factor 
in  washover  sedimentation.  With  a  storm  surge  sufficient  enough  to  keep  the 
washover  area  submerged,  flow  characteristics  will  vary;  therefore, 
textural-structural  response  properties  will  vary.  On  a  smaller  scale, 
temporal-spatial  variations  in  the  local  setting  (e.g.,  slope,  initial  flow 
velocity,  sediment  size  and  sorting,  sediment-fluid  concentration,  and 
percolation)  may  also  result  in  process-response  variation. 

The  occurrence  of  washover  suggests  an  important  relationship  to  the 
littoral  system  in  subaerial  sediment  storage  versus  permanent  offshore 
sediment  loss  or  temporary  offshore  storage.  Washover  contributes  to  the 
landward  and  vertical  accretion  of  coarser  detritus  (relative  to  other 
facies)  in  the  barrier  environment.  The  sediment  is  stored  until  later 
erosion  by  similar  storm  processes,  wind,  incremental  encroachment,  and 
erosion  by  surrounding  water  bodies  or  other  sediment  processes.  The 
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landward  sinks  are  considered  temporary  in  an  eroding  beach.  As  the  beach 
zone  encroaches  upon  the  washover  deposits  or  redistributed  washover 
sediment,  the  material  is  released  to  the  littoral  system.  (Author  Summary) 


US  Army  Engineer  Division,  North  Atlantic.  1963.  "Report  on  Operation 
Five-High,  March  1962  Storm,"  New  York,  N.Y. 

This  report  offers  one  the  most  complete  accounts  of  storm 
characteristics  and  effects  of  the  6-8  March  1962  northeaster  upon  the  US 
Atlantic  coast,  from  Virginia  to  New  York.  The  storm  lasted  through  five 
consecutive  spring  high  tides,  causing  record  water  levels,  extensive  beach 
erosion,  and  damage  to  coastal  structures.  This  report  includes  detailed 
discussion  of  storm  characteristics  and  damage,  emergency  operations,  and 
poststorm  restoration  projects.  Extensive  appendixes  include  excellent 
ground  and  aerial  photography  along  the  beaches,  beach  restoration  data,  and 
maps  of  overwash  extent  along  Long  Island,  New  York.  (Knowles) 


Wilby,  F.  B. ,  et  al.  1939.  "Inspection  of  Beaches  in  Path  of  the  Hurricane 
of  September  21,  1938,"  Shore  and  Beach,  Vol  7,  pp  43-47. 

Ground  and  air  inspection  of  some  of  the  beaches  from  Jacob  Riis  Park, 
Long  Island  to  Chatham,  Massachusetts,  was  conducted  to  record  the  effects  of 
the  hurricane  and  determine  whether  higher  and/or  wider  beaches  provide  storm 
protection.  Volume  of  sand  fills  (washovers)  appeared  to  be  greater  than 
beach  erosion  volume  in  several  localities,  possibly  indicating  transport 
from  underwater  terraces  (bars?)  or  offshore  shoals,  evidenced  by  cobbles 
with  fresh  kelp  stalks  attached  which  were  found  above  the  normal  high-water 
line.  Existing  inlets  were  widened  by  the  storm,  and  several  new  inlets  were 
formed.  A  major  geographic  effect  of  the  storm  on  Fire  Island,  New  York,  as 
viewed  from  the  air,  was  widening  of  the  island.  Many  washes  (overwash 
channels)  and  several  temporary  inlets  appeared  to  have  occurred  at  places 
where  dunes  had  been  leveled  for  construction  sites.  Generally,  the 
hurricane  had  the  effect  of  eroding  the  foreshore  and  increasing  the  width  of 
beaches  by  deposition  of  large,  wide  fills  (washovers)  landward  of  former 
berms.  Dunes  of  sufficient  height,  stabilized  by  vegetation,  provided  some 
protection  to  landward  areas.  Wider  beaches  may  have  provided  increased 
storm  protection  in  other  areas.  (Knowles) 
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